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Now for the first time, Bell Labora- 
tories’ new Optical Gas Maser con- 


Laboratories. 


The Optical Gas Maser (above) was first demonstrated at Bell Telephone 
Heart of unit is a 40-inch tube containing helium and 


tinuously generates light waves that 
are “coherent.” That is, the waves 
move in phase as seen looking across 
the beam. With further research, it 
is expected that such beams can be 
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be projected very precisely through 
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Research with coherent light is 
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beam of infrared light that may one day be used in communications. 
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ON THE MULTIDIMENSIONAL PREDICTION AND FILTERING 
PROBLEM AND THE FACTORIZATION 
OF SPECTRAL MATRICES 


BY 


E. WONG! AND J. B. THOMAS ? 


ABSTRACT 


Since the linear prediction and filtering of multiple correlated stationary random 
signals involves the solution of a matrix integral equation, several methods of solution 
of this equation have been developed. In this paper the direct solution by complex 
variable techniques is extended from the one-dimensional case. The problem reduces 
to the factorization of a spectral density matrix and this factorization is the principal 


topic of the paper. 

Under very general conditions, a factorization procedure has been developed by 
Wiener and Masani (1)’ for discrete time series. Their technique is applied here with 
some modifications to the factorization of the spectral density matrices of continuous 
processes. A simplified procedure is developed also for the case where the elements 
of the spectral density matrix are rational functions of frequency. Examples illus- 
trate the general technique and the simplified procedure. 


INTRODUCTION 


In the theory of the linear prediction and filtering of correlated 
stationary random signals, the form of the desired multidimensional 
network is specified by the solution of a matrix integral equation (2) 


(1) 


lA 
8 


where R(7r — a) is a matrix containing as elements the covariance func- 
tions of the received waveforms e;(t), that is, 


2 Associate Professor of Electrical Engineering, Princeton University, Princeton, N. J. 
3 The boldface numbers in parentheses refer to the references appended to this paper. 


(Note—The Franklin Institute is not responsible for the statements and opinions advanced by contributors in 
the JouRNAL.) 
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Rix(r7) = E}e;(t)ex(t + T) Sy (2) 
a a 1 Research staff member, IBM Research Center, Yorktown Heights, N. Y. 


88 E. anv J. B. Tuomas IJ. F. 1. 


Further, if m,(t) are the desired outputs, then the outputs of the op- 
timum network are given by 


m*(t) = f (3) 


and the elements of R,,(7) by 
= Efe;(t)m,(t + 7)}. (4) 


Thus the matrices occurring in Eq. 1 are defined by Eqs. 2, 3, and 4, 
with the notation ~ denoting transpose in Eq. 1. 

The method of undetermined coefficients was proposed originally 
by Wiener (2) for the solution of Eq. 1. An alternate method of solu- 
tion (3) is to reduce Eq. 1 to a system of simultaneous differential equa- 
tions with constant coefficients, and represents a generalization of the 
work of Zadeh and Ragazzini (4). It is the purpose of this paper to 
call attention to a third technique which may not be widely known. 


A MATRIX WIENER-HOPF APPROACH 


In a manner analogous to the one-dimensional case (2, 5) Eq. 1 can 
be rewritten as 


with 
f(r)=0, (6) 


The Fourier transforms of R(7), h(7), R,.(7) and f(7) may be denoted 
by ®(iw), H(iw), ®,,(iw) and F(iw), respectively, for example, 


Taking transforms of both sides of Eq. 5 results in 
®(iw)H (iw) — ®,,(iw) = F(iw). (7) 


It follows from Eq. 6 that F(7w) is analytic and bounded in the upper 


half w-plane (6). 
Now, let ®(zw) be factored into the form 


®(iw) = Wt (iw) W(tw), (8) 


where 
wt (iw) A w(—iw). (9) 


: 
‘ 
Racy 
( 
: 
| 
3 
ate 
: 
: 
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The factorization given by Eq. 8 is to be such that the matrix W (zw) is 
element-by-element analytic and bounded in the lower half w-plane, 
and its determinant Det W(iw) is free of zeros in the lower half w-plane. 
These two properties of W(iw) imply that 4 (zw) and its inverse W-! (zw) 
both represent transfer functions of physically realizable multidimen- 
sional networks.‘ The problem of realizing the factorization of Eq. 8 
will be deferred until the next section. 
Multiplying Eq. 7 from the left by (zw) yields 


Ww (iw)H (iw) — [wt (iw) = (10) 


Now, since both W(7w) and H(iw) are analytic and bounded in the lower 
half w-plane, the time response of W(7w)H (zw) must vanish for negative 
time, that is, 


ty 
[ W(iw)H (tw) Je"™dw = 0, 7 <0. (11) 
Similarly, it follows that 
wt (iw) F(iw) fe*"dw = 0, > do. (12) 


Therefore, the time response of Ww (iw)H (iw) must be equal to the non- 
negative-time portion of the time response of [ Wt (iw) (tw). This 
fact can be expressed by the relationship 


W (iw)H (iw) = if (iu) (iu) du dr. (13) 
TJ0 
From Eq. 13, the transposed transfer matrix H (iw) is given by 
H(iw) = f e~ier f (iu), dr, (14) 


which represents a generalization of the one-dimensional solution (5). 


THE FACTORIZATION OF THE SPECTRAL DENSITY MATRIX 


It can be seen that an explicit solution of Eq. 1 can be found pro- 
vided that the factorization (Eq. 8) of the spectral density matrix ® (iw) 
can be effected. One factorization technique has been developed re- 
cently by Youla (7). Youla’s technique is applicable whenever ®(iw) 
is element-by-element a rational function of w. An earlier factorization 
procedure was developed by Wiener and Masani (1) for discrete time 


4 Thus (iw) may be said to be the transfer function of a minimum-phase multi-dimen- 


sional network. 


: 
4 
~ 
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series. Their technique is applied here with some modifications to the 
factorization of the spectral density matrices of continuous processes. 
The procedure will now be explained and illustrated. 

The matrix (iw) to be factored is assumed to satisfy the following 
conditions : 


(1) ®(tw) is Hermitian on the real w-axis, that is, 


= (iw) 
for w real. 
(2) The determinant Det (7w) satisfies the inequality 


|In| Det ®(7w) | | 
f dw < (16) 
+ 
where | | denotes the absolute value, and k is any real finite positive 


constant (see Appendix). 
(3) Let the smallest and largest eigenvalues of ®(7@) on the real 
w-axis be \; and Xe, respectively, then 


(17) 


Conditions (1) and (2) are satisfied for all physical processes. Condi- 
tion (3) is satisfied if the spectral density matrix ®(iw) is determined 
experimentally, as shown by Wiener and Masani (1), but need not be 
satisfied by theoretically postulated spectral density matrices. 


We note that the diagonal terms ®;;(7w) of the matrix ®(zw) repre- 
sent auto-spectral densities, and can be factored easily into the form 


® = V;;(—tw)V;; (tw), (18) 


where the V;,(%@) are analytic, bounded and free of zeros in the lower 
half w-plane. Now let (zw) be a diagonal matrix with elements 


(19) 


It is apparent that the matrix ®’(7w), defined by 


(iw) = (tw) (iw) (tw), (20) 


has diagonal terms which are each equal to unity. Let the matrix 
M (tw) be obtained by subtracting a unit matrix 1 from ®’ (iw); that is, 


M (iw) ="@' (iw) — 1, (21) 


AG 

(15 

0 < 
| 
| 
o< 

| 
gh 
jk (1) A 1 
10)) jke 
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The matrix M(7w) can be written in the form 
M (iw) = My + M, (tw) + M_ (iw), (22) 


where M, is a constant matrix, M, (zw) is element-by-element analytic 
and bounded in the lower half w-plane and M_(iw) is element-by- 
element analytic and bounded in the upper half w-plane. It is recog- 
nized that these requirements do not uniquely determine the decom- 
position (18). However, it is shown in the Appendix that different 
decompositions satisfying the analytic requirements yield equivalent 
results. In practice it is frequently convenient to add the require- 
ment that 

M,(«) = M_(o) = 0, (23) 


which makes the decomposition (18) unique. Now let the matrix 
N (iw) be defined by the series 


N(iwv) 41 —M, + (MM,), —[M(MM,),], 


where the subscript + has the same meaning as in Eq. 22. It follows 
immediately from the work of Wiener and Masani (1) that 


(1) The matrices and N-'(tw) are element-by-element an- 
alytic and bounded in the lower half w-plane. 
(2) The matrix G, defined by 


G A N' (iw) (iw) N (iw), (25) 


is a constant symmetric matrix. 
Furthermore, it is seen from Eq. 24 that 


Nt(o) = N(o) = 1, (26) 

if Eq. 23 applies. Therefore, 
G = @'(~) (27) 
since G is a constant matrix. Thus G can be found without performing 


the matrix multiplications implied by Eq. 25. 
From Eqs. 20 and 25, ®(iw) can be written as 


®(iw) = Nt (iw) (iw) (iw). (28) 


All that remains to complete the factorization is to factor G. This 
factorization is easily accomplished since G is a constant matrix. For 
example, G may be diagonalized by unitary operations and the square 
root taken of the resultant diagonal matrix. However, there is no need 


: 
: 
fy 
‘ 
: 
+ 
A 
: 
: 
ig 
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to proceed further than Eq. 28. Using arguments similar to those 
leading to Eq. 14, it is easy to show that, instead of Eq. 14, the solution 
can be expressed as 


~ 


SiGe) = dre-* 


x | (iu) (29) 


A SIMPLIFIED PROCEDURE FOR RATIONAL SPECTRA 


Since Eq. 24 is an infinite series, N (7m) cannot always be obtained in 
closed form. If ®’(zw) is element-by-element a rational function of zw, 
an alternate procedure yields N(zw) in a finite number of steps. 

To derive this procedure, ®’(7) is written as 


®'’ (iw) = 1+ My, + M, (tw) + M_(tw). (30) 


It should be noted that the successive terms in the series (24) for N(¢w) 
can be written more explicitly ; for example, 


(MM ,), = M.M, + M,M, + (M_M,).. (31) 


If the product ®’ (tw) N (zw) is formed and use is made of Eqs. 30 and 31, 
it is found that most terms cancel and the product becomes 


®' (iw)N (iw) = 1 + M, + [M_(iw)N (iw). (32) 


From Eq. 32, N(zw) can be written formally as 


N(iw) = [@’ (iw) + My + [M_(tw)N (iw) ]_}. (33) 


Since N(iw) is analytic and bounded in the lower half w-plane and 
[M_(iw)N (iw) |_ is analytic and bounded in the upper half w-plane, the 
poles of N (zw) must result from the upper half w-plane poles of [@’ (zw) }-'. 


Let [’(t#) }-' and M_(iw) be written as 


[ (tw) = + + tw) 


m= 


: 

: 
: 
! 

as 

Ste 

: 

. 

4 

+ > A, (3 
: and 
M_(iw) = B,.(7, — 35 
: 
3 n=1 

: 
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From Eq. 34 it is apparent that N(zw) can be written as 


K 
N(iw) = 1+ + (36) 
m=1 
Substituting Eqs. 34 through 36 into Eq. 33 and equating coefficients 
of the upper half w-plane poles yields 


L 
= A,,(1 + M,) + (on, + Tn) 
=1 


L K 
4- A. B,C, 4. -+- Tn), m K. (37) 


n=1 k=1 


Equation 37 represents a set of K matrix equations from which the K 
unknown constant matrices C,, can be found. 


Example I 

A two-dimensional prediction problem is considered in some detail. 
Let two correlated signals e;(7) and e2(7) be received in the interval 
—« <7r<t. It is desired to estimate e,(f +a) and e.(¢ + a) for 
a >0. The spectral density matrix of e:(7) and e2(7) is given as® 


1 
1+ (1 — tw)? 
lel 
€ 1 | 


(1+)? 1+ 
The matrix ®,,(7w) for this case is simply 
®,, (iw) = e**@(iw). 


From Eq. 20 it is found that 


(40) 

1+ tw 
In this example the power spectral density matrix is a rational 
function of zw; hence the simplified procedure previously described can 


be applied. For comparison the problem is solved using first the general 
approach and then this simplified procedure. 


5 This example appeared in (2) and was solved by the method of undetermined coefficients. 


9 3 
; 
4 
®(iw) = (38) 
(39) : 
1 
1 ; 
4 
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The unit matrix 1 may be subtracted from Eq. 40 to yield 


1+ iw 
Q € 
| 1 — Ww 
M (iw) = | (41) 
1 0 | 
(| 1+ 7 
and, finally, 
Q) O | 
M.(iv) =| 5, (42) 
1+ 


For this simple problem the series (24) for N(w) terminates after two 
terms, and the result is 


N(iw)=| |. (43) 
— | 
1 + tw 


From Eq. 25 the matrix G is given by 


4 
1 — iw | | 1 — tw | | 
G = | | | | 
| | 
1+ tw L1+ 
1 
G = (44) 
1 


It is seen that G is indeed a constant matrix. From Eq. 29 and 39 
the solution is found to be 


H (iw) (tw) N (tw) are | (45) 
0 LT 


Following the simplified procedure of the previous section, Eqs. 30 
and 40 are used to obtain 


5 

: 

1 Q | 

a 

> 
1 0 

2€a 

4 

: 
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(x) | M, = 
The matrix given by Eq. 40 may be inverted to yield 


1+ Ww 
€ 


1 — Ww 


1 — w 
+ iw 


A comparison of Eqs. 34 and 47 shows that a; = 1 and that 


From Eq. 41, M_ (iw) is given by 
M_(iw) = 


and Eq. 35 yields 7, = 1 and 


(50) 
0 


Substitution of Eqs. 46, 48, and 50 into Eq. 37 yields a matrix equation 
which may be solved for C, to give 


(51) 


From Eq. 36, N(7w) may be written as 


0 O 
N(iw) = 1+ 


which is identical to Eq. 43 obtained by the general procedure. 


: 
i 95 
(46) : 
(iw) = —— (47) 
l-—e 
| 
0 0 } 
— Ze 
2e 
(49) 
0) 0 
0 | 
0 
(52) 
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Example II 


In this example a somewhat more complicated spectral density 
matrix will be factored. Let this matrix @(7w) be 


9 


ot €aB 
w (a@ + (6 — lw) 
| 
| 


®(iw) = | < 1. (53) 
B° 


(a — tw)(B + iw) B? + 


The matrix (zw) of Eq. 19 is found, by factoring the diagonal terms, 
to be 


0 | 
| 
(iw) = | (54) 
B+ tw | 
B 


From Eqs. 20 and 54, @’(zw) is found to be 


a — + tw ) 


1 € | 
— tw | 
®' (iw) = (55) 
at+iwB — Ww 
— wBt+ tw 


Thus M (iw) of Eq. 21 is given by 


a — wB + tw ) 


€ 
a + iw B — tw | 


M(w) = 


| 
| 
a+iwB — tw 


€ 
lL a—twB+tmw 


From Egs. 22, 23, and 56 the matrix M, (zw) is 


2a B-a 


0 €- 
at+mPpBp+a 


M.(iw) = 
28 
+ twa 


The series (24) for N(zw) in this example is a geometric series, and can 
be evaluated to yield 


: 

: 

Z 

) 

7 
f 

- 

(56) 

omy 

(9/ ) 

0 

2 

ii 
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ky 
a + tw 


ie! 


where the constants k; and k» are given by 


= — 8) 
+ B)? — 4€ap’ 
and 
2e(a? — 
(a + B)? — 


Application of Eq. 27 yields immediately 
(61) 


The matrix G could be obtained also from Eq. 25. Since this example 
has a rational spectral density, it could be solved by the simplified pro- 
cedure previously discussed. 


APPENDIX 


In order that continuous processes may be brought within the framework of the work of 
Wiener and Masani (1), a mapping of the w-plane into some z-plane must be effected such that 
the real axis of the w-plane is transformed into the unit circle in the z-plane, the upper half 
w-plane into the region outside the unit circle in the z-plane and the lower half w-plane into the 
region inside the unit circle. Such a mapping can be obtained by the transformation 

so O<b< (62) 
Physically, the constant & can be interpreted as a scale factor which translates the dimension- 
less time scale of the discrete case into the time scale of the continuous case. It will be shown 
that the transformation (62) also serves to give a precise definition to the decomposition (22). 
Let F(z) be a function which can be represented by the expansion 


(63) 


n 


Then the functions F,(z), F_(s), and Fo are defined as 


Fy (64) 


F(2)4 3 (66) 


n=l 


; 
Satan : 
Qa 
2 
a + tw 
bie N (iw) = (58 
’ 
(59) 
(60 
i 
Le 1 
| 
2 
é 
on 
= 
Cape ; 
: 
3 
: 
= 
24 
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It is apparent from Eqs. 65 and 66 that 


F,(0) = =0. (67) 

From Eq. 62 it is seen that the points z = 0, © correspond to w = — tk, ik, respectively. 
Therefore, the condition to be imposed on the decomposition (22) is 

M,(—ik) = M_(ik) = 0, 0<k< (68) 


It may appear from Eq. 68 that the condition (23) is inapplicable, since k is restricted to 
be finite. This, however, is not the case. Consider the alternative decompositions of M(iw) 


M(iw) = M,(iw) + M_(iw) + Mp = M* (iw) + M~(iw) + (69) 


where M, and M* are analytic in the lower half w-plane, M_. and M~ analytic in the upper half 
w-plane and My and M’ are constants. Furthermore, require that 


M,(—ik) = M_(ik) = 0, 


0<k< ~, (70) 


and 


Mt(o) = M-(0) = 0. (71) 


It is apparent that M* and M,, thus defined, differ only by an additive matrix constant. Now 
let N(iw) be given by 


N(iw) = 1 — M, + (MM,), — [M(MM,),], + --- 


and let N’(iw) be given by 


N’(io) = 1 — M* + (MM*)* — [M(MM*)*}* + 


Equation 73 can be rewritten as 
N’(iw) = 1 — (M, + K,) + {[M(M, + K,)], + K.} --- 

{1 — M, + (MM,), — [M(MM,),], + ---}{1 —K, + K. —K;+---} (74) 
N(iw)K. 


Thus the matrices N’(iw) and N(iw) differ only by a multiplicative constant matrix. This 
difference does not affect the final result given by Eq. 29. To see that such is the case, instead 
of Eq. 29, write 


From (74), it is found that 


(N’(iw) = K'N' (iw), (76) 
G’ = [N'(iw) ®'(iw) N’'(iw) 
= (iw) @®'(iw)N(w)K 
= K'GK, (77) 
and 
[G’}? = (78) 


Therefore, 


x J, {5 K (iu) iu) | (79) 


> 


{1 


) 
2a 


dre-‘r 


= H(iw). 


ae 
4 
: 
: 
(72) 
(73) 
\ 
: 
: 
: 
f 
4 
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SOME NEW CONCEPTS IN NONLINEAR SYSTEMS 


BY 


C. LAKSHMI-BAI! 


ABSTRACT 


Some new concepts characterizing the response of nonlinear systems are de- 
veloped. These new concepts are denoted by the terms, the transient system 
equivalent, the response vector, and the space-phase components. This third concept 
is analyzed in comparison with the well-known technique of symmetrical components. 
The performance of a multiplicative feedback control system is represented by a 


nonlinear integro-differential equation; its solution is obtained by the principle of 
variation of parameters. The system response is treated as a vector and is resolved 


into its space-phase components. The individual effects of these components on the 


performance of the system are discussed. The suitability of the technique for the 


transient analysis of higher order nonlinear control systems is discussed. 


INTRODUCTION 


Of late there has been a lot of active interest displayed in the general 
field of nonlinear systems. Various techniques have been developed to 
analyze different aspects of nonlinear behavior. The phase-space 
technique originally developed by Cauchy (1)? and subsequently ap- 
plied by Ku (2) to physical problems is an elegant method for the 
purpose of analysis. The describing function technique is a useful 
design tool for analyzing the system behavior under steady state. 
The assumptions necessary and the accuracy attainable by this method 
are discussed in detail by Johnson (3). Chen (4) used the principle of 
quasi-linearization to arrive at the transient response of the systems 
when the nonlinear components are of the piecewise type. Under 
many restrictions, a similar method has been developed by Grensted (5) 
for the transient analysis of second order nonlinear control systems. 
The purpose and aim of this paper are to advance a technique which 
offers a new perspective and allows greater insight into the physical 
significance of higher order nonlinear systems. In the development of 
this method, based essentially on the principle of the variation of 
parameters, some new concepts have been recognized. These concepts, 
when investigated in detail, led to very interesting properties. The 
most important feature of the concept of space-phase components is its 
striking analogy to the classic technique of symmetrical components (6). 
Details of working in case of a third order nonlinear control system 
are illustrated, and the results are compared with those obtained by 


1Senior Research Assistant, Department of Power Engineering, Indian Institute of 
Science, Bangalore, India. 
2 The boldface numbers in parentheses refer to the references appended to this paper. 
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computer analysis. The concepts are extended to fourth and higher 


order systems. 


NOMENCLATURE 


r(t) = Reference input to system 
e(t) = instantaneous error 

C(t) = system response 

y(t) = instantaneous phase of oscillation 
w(t) = instantaneous frequency 


II 


dy 
— time rate of change of phase, 


dt 
multiplying unit 


BASIC THEORY 


The solution of a nonlinear differential equation of a given order is 
closely related to the number and the nature of roots of the associated 
linear differential equation. Further, the solution of an mt order 
differential equation contains m arbitrary parameters which can be de- 
termined by utilizing the 2 boundary conditions available. 


Definitions of New Concepts 


Let the response vector be defined as the solution of the given 
nonlinear differential equation. The significance of the term ‘‘re- 
sponse vector’ is to indicate the interconnection of a physical behavior 
with a mathematical analysis. The term “‘response”’ itself stands for 
the response of the physical system, while the term ‘“‘vector’’ stands for 
the corresponding mathematical analog. Thus the degrees of freedom 
of the physical system represented by the order of the differential equa- 
tion corresponds to the number of components of the vector, or the 
dimensions of the space to which the vector belongs. In other words, 
while the dimensional space of the response vector depends on the order 
of the governing equation, the nature of its components varies with the 
type of input and the associated initial conditions. 

Let the ‘‘transient system equivalent”’ be defined as a mathematical 
representation of the original nonlinear system. In such a representa- 
tion there are no more cross-harmonic terms present, yet the original 
nonlinear effects have been adequately accounted for by their d.c. and 
first harmonic terms. This mathematical representation takes the 
form of a vector referred to an orthogonal reference frame in a space of 
as many dimensions as the order of the original differential equation. 
As will be seen later, this representation is valid in the transient as well 
as in the steady state. . 

The two concepts as defined above naturally lead to the third con- 
cept of “‘space-phase components.’’ These are the projections of the 
response vector along the various co-ordinate axes of the orthogonal 
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reference frame. A geometrical interpretation is possible in case of a 
third order system, and the concepts can be extended to the analysis 
of higher order systems. This to some extent explains the need for the 
terms ‘‘space’’ and ‘“‘components.’”’ The term ‘phase’ also has its 
physical significance, based on the choice of the reference axes, along 
the directions of the in-phase, the quadrature and the zero-phase com- 
ponents of the response vector. The in-phase is with reference to the 
direction, predetermined on the basis of the associated linear differential 
equation, while the quadrature component is at right angles to the 
in-phase component. The zero-phase component, on the other hand, 
includes all the terms which are independent of phase. Such terms may 
arise from the real roots of the associated linear differential equation, 
as well as the d-c. components due to the additional nonlinear terms. 
Thus, it can be easily visualized that the direction of the zero-phase 
component is at right angles to the plane containing the in-phase and 
the quadrature components. Further, the order of the zero-phase 
component is equal to the number of real roots of the associated linear 
differential equation. For example, in case of a third order system, the 
order of the zero-phase component can be one, or three as the case may 
be, while in the case of a fourth order system, it can be two or four, 
because complex conjugate roots have to appear in pairs. 

Thus, one can correlate the degrees of freedom of the response vector 
with the order of the governing differential equation. Further, the 
nature of the space-phase components is tied up with the nature of the 
roots of the associated linear differential equation, modified to some 
extent by the nature of the non-linear terms present. 

For the sake of simplicity, the following discussion is with reference 
to a third order system, though similar principles apply to higher order 
systems, with a few additional remarks, which are enumerated in 
later pages. 

For a moment let one consider the total degrees of freedom of three 
vectors in a three-dimensional space. With no restrictions whatever, 
each of the vectors has three degrees of freedom; thus the three vectors 
together have nine degrees of freedom. However, let the following two 
restrictions be imposed. One of the three vectors is along a specified 
direction (determined by the phase of the assumed solution) and the 
other two directions are mutually orthogonal to it. Now each of such 
vectors has only one degree of freedom, which is influenced by the 
corresponding magnitude. Thus, three such vectors together have 
only three degrees of freedom. This statement can be easily verified 
with the aid of spherical coordinates. Hence, it is possible to represent 
the solution of a third order differential equation as a response 
vector in a three-dimensional space, having three mutually orthogonal 
components. 
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Analogy with Symmetrical Components 


It is well known that by the method of symmetrical components, a 
set of unbalanced currents (or voltages) may be resolved into systems 
of balanced currents (or voltages), equal in number to the number of 
phases involved. The principle in essence is that an unbalanced multi- 
phase electrical system can be represented by a simple mathematical 
analog. Among the varied advantages of this technique, those that are 
of interest to the present work are the stability study of power systems, 
determination of responses under unbalanced fault conditions, and 
evaluation of time constants. Table I shows the principal analogs 
between symmetrical components and the space-phase components, as 
developed in this paper. Table II shows some aspects of symmetrical 


TABLE I.—Analogs. 


SpAcE-PHASE COMPONENTS 


SYMMETRICAL COMPONENTS 


Number of phases of the electrical system Order of the governing equation 


Nonlinear effects 


Unbalance effects 


Mutual coupling between phases due to — Presence of cross-harmonic terms 
dis-symmetry 


Positive and negative sequence com- — In-phase and quadrature components 


ponents 


Zero-sequence component Zero-phase component 


A simple mathematical representation of | A simple mathematical representation of a 


an unbalanced electrical system nonlinear physical system 


~ 


Enables stability study of power systems Enables stability study of physical systems 
governed by the nonlinear differential 


equation 


oo 


Valid for analysis of steady state, in its Valid for analysis of transient as well as 
original form. Transient analysis can steady state behavior 


be effected by the modified Lyon's com- 


ponents 


Various time constants of the systems can ‘Time of rise, time of settling and other rele- 


be determined vant time constants can be determined 


Initial currents and voltages can be ac- Initial conditions can be accounted for 


counted for 


Suitable for the study of unbalanced fault Any type of input in closed form, can be 


conditions studied, and instruments can be devised to 


control extraneous disturbances 


Grounding phenomena associated with Zero-phase component influences only the 

zero-sequence component amplitude of the response, while the 
in-phase and quadrature components in- 
fluence both the amplitude and_ the 
frequency 
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h 
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TABLE II.—Differences. 
SYMMETRICAL COMPONENTS Space-PHASE COMPONENTS 


Linear properties are assumed in the deriva- | The only assumption made is that higher 
tion of the equivalent representation harmonics in the response of the nonlinear 
component are adequately filtered out by 
the linear part. This assumption is quite 
valid and practical with nonlinear control 

systems 


Sequence currents, voltages and impedances ‘The only variable is the system response 
Evaluation of sub-transient parameters None 


Synthesis of networks to represent faults Synthesis of components to represent the true 
nonlinear effects is not possible because of 
the assumption made 


5 Special features in the application to rotat- No essential difference in the application to 
ing machines static and dynamic nonlinear components 


Sequence components are coplanar for any Space-phase components are mutually orthog- 
number of phases onal and are located in space 


components for which there is no counterpart in the space-phase 
components. 


Fic. 1. Multiplicative feedback control system. 


Contrast with the Phase-Space Technique 


The phase-space technique portrays the relation between the system- 
response variable and its derivatives as a surface in a space of dimen- 
sions equal to the order of the governing equation. Though the space- 
phase components are also equal in number to the order of the governing 
differential equation, this technique represents the resolution of the 
response vector into mutually orthogonal components. The phase- 
space technique is topological in nature while the method of this paper 
is analytical. The phase-space technique is absolutely rigorous, while, 
the method of this paper is valid under one assumption (see Table IT, 
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item 1), which is practically justified with closed-loop nonlinear control 
systems. Thus, it can be concluded that the two methods are entirely 
different in their principle and application. 


Illustrative Example 


In Fig. 1 is shown a multiplicative feedback control system which 
has been studied by Ku and Chen (7). In Fig. 2 is shown the equiva- 


Fic. 2. Equivalent single-loop system. 


lent single-loop system. The governing differential equation of this 
system, to unit step input is given by 


é+é+(i1- be =: (1) 


where e stands for the instantaneous error and the dots denote differen- 
tiation with respect to time. The initial conditions which the system 
has to satisfy are given by 


e(0O) = 0.5, é(0) = 0, and é(0) = 0. (2) 


Since the principal aim of the present analysis is the application of 
some new concepts, only the stable case as determined by Ku and Chen 
is taken up for investigation. Further, transient analysis of systems is 
of interest only when the system is known to be stable, as determined 
by frequency response (8) or other preliminary design criteria. The 
associated linear differential equation is found to have one real root 
and a pair of complex conjugate roots. Therefore, one might assume, 
by using some engineering judgment, that the solution of Eq. 1 could 
be oscillatory in nature. Thus, a solution of the type shown in Eq. 3 
may be assumed, 


e(t) = a(t) sin Y(t) + d(?). (3) 


Equation 3 indicates that the frequency of oscillation is itself time- 
varying. Such, in fact, is a special feature of nonlinear systems. The 
actual solution in the forrn of Eq. 3 will be based on the principle of 
variation of parameters. 
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As a first approximation a is assumed to be a constant equal to do, 
so that the solution becomes, 


(4) 


ei(t) = ay sin y(t) + dit). 


Further assumptions made in the first approximation are : 


1. Since e(t) represents the instantaneous error, in a stable system, 
the amplitude of the error is decreasing with time. Hence one may 
assume |a| < 1, so that a? and higher powers of a may be neglected. 

2. The frequency is slowly varying over the length of one cycle, and 
therefore one may neglect the second and higher derivatives of y, 
the phase. 

3. b(t), due to the real root, has a strength less than 1 and is also 
slowly varying. Thus higher powers of }b, and its second and higher 
derivatives may be neglected. 

4. The higher harmonics in the solution for e,.or the response vector, 
may also be neglected, since the linear part will act as an effective low- 
pass filter. 

5. Higher powers of ¥, the frequency, may not be neglected, as it 
is likely to be greater than 1. 


It is important to observe that the precise basis of the above assump- 
tions is not always too clear. Their validity can only be confirmed 
after comparison with computer analysis. However, an experienced 
designer is capable of exercising this discretion, and the likely errors 
are expected to be within tolerable limits. 

Substituting from Eq. 4 into Eq. 1, and simplifying according to the 
above assumptions, one gets 


i 
(— ao¥? + — aobb) sin + (any — av¥*) cosy + (4b + b) = 0. (5) 


As discussed under the basic theory, coefficients of each of the three 
phase components in Eq. 5 can be set equal to zero, thus leading to 
three simultaneous equations: 


doy = 0 
b+ =0. 


Solution of Eqs. 7 and 8, followed by further integrations leads to 


b(t) = 


y(t) = ¢t + constant. (10) 


on 
Me 
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(8) 
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Because of the first approximation, Eq. 6 is of little use. 
Using the initial conditions according to Eq. 2, Eq. 10 reduces to 


=t+>. (11) 


Now, substituting in Eq. 4 from Eqs. 9 and 11, one gets 


= do sin (: + =) + (12) 


As a second approximation, let a be a function of time a(t), so that 


e,(t) = a(t) sin (: 4 + (13) 


By physical considerations a(t) may be assumed to be exponential in 


nature. 
Thus, let 


a(t) = doe. (14) 
Using Eqs. 13 and 14 along with Fas. 1 and 2, one gets, 
e(0) = ao + by = 0.5 (15) 
é(0) = aay — 0.5d5 = 0 (16) 
E(O) = — ao + = 0. 
Solving Eqs. 15, 16 and 17 simultaneously, one gets 
ao = — }, and by = 2. 


Now, the required solution according to Eq. 3 is given by 
e(t) = — ge*'sin > + (19) 


Checking for the accuracy of the assumptions made with the aid of the 
specified initial conditions, one finds that, é(0) = — 4, instead of zero. 
For a unit step input, Eq. 19 leads to the system response, 


C(t) = 1 + sin (: + — (20) 
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As shown in Fig. 3, the evaluated analytical solution for Eq. 20, 
agrees very favorably with that recorded on the computer (2, p. 320, 
Figs. 12-24). 

In the above derivation, the solution for e(¢) assumed on the basis 
of the associated linear differential equation is of the type 


e(t) = sin W(t) + doe®*. (21) 


c(t) 


“4 5 6 7 8 
TIME IN SECONDS. —> 


Response curve from analog study 
*oints calculated from analytical study. 


Fic. 3. 


Such a solution contains five parameters to be determined from the 
three initial conditions according to Eq. 2. The additional equations 
necessary are supplied by utilizing the concept of the space-phase 
components. 

Thus the final solution is obtained by the process of successive 
approximations based on the principle of variation of parameters. 

Comparing Eqs. 21 and 20, one finds that the zero-phase component 
of the response vector is given by 1 — 3¢-“?, while the in-phase com- 


ponent is given by §e-*! sin (: +z ) ; the quadrature component itself 


is zero. It can be easily visualized that the zero-phase component is 
of order one and modifies only the amplitude, while the in-phase com- 
ponent effects both the amplitude and the frequency. 


DISCUSSION 


During the course of the solution, the coefficients of terms contain- 
ing sin y, cos y and terms independent of ¥ are separately set equal to 


0.54 
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zero. Krom the point of view of equilibrium of forces, this stands for 
the fact that the components along mutually orthogonal directions of 
the resultant force are individually zero. This further led to the recog- 
nition of the fact that the net response vector can be resolved into space 
components which characterize the phase properties. It has also been 
observed that the response vector has already accounted for the non- 
linear effects, and therefore such a resolution into space-phase com- 
ponents is independent of the nature of the system. 

The technique developed in this paper is particularly suitable for 
multiplicative feedback control systems, because asymmetric and 
dynamic nonlinearities do not present any special difficulties. The 
additional advantages of the method as compared to the describing 
function technique are that different types of inputs in closed form can 
be effectively treated, along with specified initial conditions. This latter 
fact makes the technique particularly useful for evaluating the im- 
portant transient characteristics such as the peak overshoot, time of 
rise, time of settling and also the possible limit-cycle behavior. 

It is to be observed that with systems higher than the third order, 
prior engineering judgment is necessary to make an approximate guess 
about the several modes of oscillation and the possible number of real 
roots. For example, in the case of a fifth order system, there may be 
only one pair of complex-conjugate roots and three real roots, which 
means that the corresponding zero-phase component is of order three. 
It may also happen that the fifth order system may exhibit two modes 
of oscillation due to two pairs of complex-conjugate roots and the 
remaining one real root contributes to a zero-phase component of order 
one. Thus the concepts discussed with respect to third order systems 
can be carried over to higher order systems in analysis. However, 
geometric representation is no longer possible for systems of order 
greater than the third. 

It is easily visualized from the above discussion, that the zero-phase 
component influences only the amplitude of the response, while the in- 
phase and quadrature components affect the amplitude as well as the 
frequency of the response. Thus it is seen that these new concepts 
isolate important characteristics into quantities which represent a better 
criterion of the controlling factors. This feature makes it a powerful 
tool in the synthesis of nonlinear systems. 


Advantages of the New Approach 


1. An approximate solution for the nonlinear differential equation 
can be assumed, based on the associated linear differential equation. 

2. The above fact enables the use of the method of variation of 
parameters. 
3. Successive approximations improve the degree of accuracy. 
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4. No special difficulties are encountered with dynamic nonlinear 
components. 

5. Particularly suitable for multiplicative feedback control systems. 

6. Mathematical representation is free from cross-harmonic terms, 
yet the nonlinear behavior has been adequately accounted for. 

7. D-c. terms arising from nonlinear components form part of the 
zero-phase component ; thus they can be easily taken care of. 

8. Resolution into space-phase components provides the additional 
equations necessary for evaluating the parameters in the assumed 
solution. 

9. Effects on the system due to the different phase components are 
isolated. This is a desirable feature in synthesis. 

10. Suitable for transient analysis of higher order nonlinear control 
systems. 


CONCLUSIONS 


1. Three new concepts—the transient system equivalent, the re- 
sponse vector, and the space-phase components—are developed. 

2. These new concepts are compared with the classic technique of 
symmetrical components. 

3. An illustrative example is worked out, discussing the various 
aspects and advantages of this new method of solution of nonlinear 
systems. 

4. It is believed that this technique will lead to further developments 
in the theory and application to nonlinear control systems, possibly 
even for random inputs. 
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RESTRICTIONS ON THE NATURAL FREQUENCIES 
OF AN RC-RL NETWORK * 


BY 


D. A. CALAHAN ! 


ABSTRACT 


The following problem is investigated: given any polynomial N(s) of the form 


N(s) = Il (s + pi) Il (s + s;)(s + &) 
1 


what are the necessary and sufficient conditions that N(s) be formed in the following 


manner: 


N(s) = AoII(s + ai) + BoII(s + Ao, Bo > 


where the —a; and —0; alternate at most in pairs along the negative real axis. Such 
formation is necessary for the —p; and —s, to be the natural frequencies of a network 
consisting of RC and RL subnetworks joined at two terminals. 

It is shown that N(s) may be so formed if and only if 


m 
(a) torn = 0, Darg s; < 
1 


(b) forn > 0, > Oand Darg s; < 2/2 
1 


Also, a simple method is presented for finding suitable a; and 6; whenever formation 
is possible. 

The result is expected to have application in the synthesis of both passive and 
active networks. 


It is well known that if a network function has complex poles, then 
any passive realization must contain at least one inductor. DeClaris 
(1)? and Hakimi (2) have studied the natural frequencies of this “single 
inductor” case (Fig. 1). Horowitz (3) has similarly obtained restric- 
tions for the ‘‘two-inductor’’ model of Fig. 2. It is now desired to 
obtain necessary and sufficient conditions that a given set of natural 
frequencies may be those of a network of the form of Fig. 3, where the 
RC and RL networks are of any complexity. 


R Ef 2 RC RL 


> 


Fic. 1. Single inductor Fic. 2. Two inductor model. Fic. 3. RC-RL model. 
model. 


* The work reported here was supported by the Air Force Office of Scientific Research 
under Contract AF(638)-63. 
1 Department of Electrical Engineering, University of Illinois, Urbana, IIL. 
* The boldface numbers in parentheses refer to the references appended to this paper. 
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THE APPROACH TO THE PROBLEM 


Let a polynomial N(s) contain all the given natural frequencies. 
Then, for these frequencies to be those of the RC-RL model, NV(s) must 
be formed in the following manner: 


= + (1) 


where Zrz(s) and grc(s) are RC and RL driving point impedances,’ 
respectively, and where D(s) is unspecified, it being composed of the 
poles of zez(s) and Zgrc(s). The alternation of the poles and zeros of 
Zrx(s) and Zrc(s) along the real axis imposes restrictions on the zeros 
of N(s). The nature of these restrictions is the subject of this 
investigation. 

To begin, Eq. 1 may be re-written 


N(s) _ mrx(s) mre(s) (2) 


D(s) drx(s) dre(s) 
N(s) nru(s)dre(s) + 


D(s) dri(s)drc(s) 
N(s) A ovAn(s) + BowBy (s) (4) 
D(s) dri(s)dre(s) 
Hence, 
N(s) = AonAn (s) > BonBy(s) (5) 
N(s) = Aon II (s + a:) + Bon II (s + b;) (6) 
ai S 
bs < 
ai <b, 


A necessary condition, then, that the zeros of N(s) be formed by the 
sum of an RC and an RL impedance is that N(s) be decomposable into 
an AovAy(s) and a BoyBy(s) whose zeros can be manipulated to form 
the component impedances. The following theorem will place more 
meaningful bounds on A y(s) and By(s). 


Theorem I: Let the number of zeros of Ay(s) and By(s) to the right 
of any point ¢ on the real axis be Z4(¢) and Z x(a), respectively. Then 


2 >Za(c) > 0 (7) 


3 The formation of N(s) from RC and RL admittances is an identical problem. 
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Proof: Let the number of zeros of mp1(s), Mec(s), dei(s), and 
drc(s) to the right of any o be Z,,(¢), Zn,(0), Za,(o), and re- 
spectively. From the properties of RC and RL impedances, 


1 > Zn,(¢) — Za,(9) 


IV 


0 (8) 


V 


Adding Eqs. 8 and 9, it is found 


2 + Za.(0) — — 2 0. (10) 


But, from defining Eqs. 3 and 4, 


Zni(a) + Za (9) = Za(o) (11) 


(12) 


and 


+ Zn(o) Z 


From Eqs. 10, 11 and 12, then, 


2 > Za(c) —Za(c) > 0. (13) 


THE OPTIMUM DECOMPOSITION 


The useful necessary conditions imposed upon AoyAwy(s) and 
BonBy(s) by the properties of the RC and RL impedances are now 
summarized : 


(14) 


. Aoy and Boy are real and positive 


(15) 


. a; and 6; are real 


3. > Za(o) Z = 0 (16) 


4. —a; and —); are non-positive. (17) 


The question of whether a given N(s) can be decomposed into an 
AonAn(s) and a BoyBy(s) which fulfill Conditions 14-17 will now be 
considered. 

The satisfaction of Condition 17 is clearly guaranteed from the 
non-positiveness of —a,. If it were possible to find a decomposition 
of N(s) which, of the class of all decompositions which satisfy Condi- 
tions 14-16, minimizes —a,, then only the —a, of this decomposition 
would have to be examined to predict the possibility of RC-RL de- 
composition of N(s). Such a decomposition will be termed the opti- 
mum decomposition and will be written N(s) = Ai(s) + BoiBi(s). Its 
form is now deduced. 
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Theorem II; Let N(s) be a polynomial with real zeros and pairs of 
complex zeros, so that 


N(s) = JI (s + II (s + 5:)(s + (18) 
1 1 
Then of all the decompositions 
N(s) => AonA n(S) + BonBy(s) (19) 


where A y(s) and By(s) fulfill Conditions 14-16, the decomposition that 
minimizes —a, is of the form 


(20) 
and 
ByBy(s) = Bor IL (s + ps) TL (s + (21) 


Proof: Givenany A ovA y(s) and BoyBy(s) which satisfy Theorem 
I, it will be shown that they can always be converted to A,(s) and 
Bo, B,(s) while —a, is decreased or remains the same. Let 


(s +a) (22) 
and 
BoxBw(s) = Bow I (s +b) (23) 
where 
2>q-r2>0 (24) 


in accordance with Theorem I... Then form 


A ony1Awyi(S) = AonAn(s) + 50 II (s + (25) 
1 
Bon 41By41(s) BonBy/(s) — bo Il (s + 5;) (26) 
1 
where 
(27) 


so that Boy 1 remains positive for some 6) > 0, and where 6; are located 
at any double zeros of Ay(s) and By(s). 
As 6» increases from zero, one of three things can happen which may 
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Zeros of A,(s), B,(s) 


$2 6) 


Production of Double Zeros of Ay, oS) and By, 


& 


Zeros of Ay (s), By (s) 


Real Zero 


Production of Real Zero of N(s) 


A X 
Zeros of Ay(s) and By (s) 


-5 3 
©- 


X X 


4 
Production of Zero at Infinity 


A. 


~— 
Zeros of Ay (s) and B,(s) 


A 


First Step 


Second Step 


(d) Steps in Reduction to Optimum 


Fic. 4. 
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cause the resulting Ay4+:(s) and By4:(s) to violate Theorem I. They 
are listed below, together with a procedure to avoid the violation. 


1. For some 69,, any zeros of Ayw4:(s) or By4:(s) could become 
double before becoming complex (Fig. 4a). If this occurs, choose 
5) = bo,. Then by use of Eqs. 25 and 26, form a new decomposition 
from Ayyi(s) and By,i(s). Note that the 6; have now changed from 
the original. 

2. For some 6, real zeros of Ay4i1(s) could become equal to real 
zeros Of By,.(s) (Fig. 40). If this occurs, choose 6) = 59, and repeat 
Eqs. 25 and 26 for 


II (s + II (s + pi) 


where —/; are the common real zeros of Aw,i(s) and By,,(s). These 
are, of course, real zeros of N(s) also. 

3. For some 6o,, if ¢ = r, a zero of By,:(s) could travel to — © 
(Fig. 4c). In this case, choose 6) = 69, and repeat Eqs. 25 and 26 for 
the resultant A and By,,(s). 

In any of steps 1, 2, or 3 it is clear that no zeros of Awy4:(s) can 
occur to the right of —a, of Ay(s) since this region is void of any zeros 
of 5(s). Therefore, —a; must be decreased or, at worst, remain the 
same. The latter will happen only when a double zero of Ay(s) occurs 
at —a,. Thenext step must then decrease —d,, as Fig. 4d demonstrates. 

For any decomposition to minimize —a,, it must be immune to 
improvement by all of the preceding steps. But the only assumption 
made in forming new decompositions is ¢ < r, which therefore must not 
be satisfied by the optimum decomposition. 

Now ¢ > +r requires that there be more pairs of double zeros of 
Ay(s) and By(s) than there are zeros of By(s). This does not include 
the real zeros of N(s) that have been omitted from the process by 
Step 2. Let 


= + (29) 
+ (30) 


r 


where g; and 7; are the number of pairs of double zeros of Ay(s) and 
By(s) and where g2 and 72 are the number of their simple zeros, re- 
spectively. Then ¢ > r demands 


(31) 
or 


} 
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Also, from Eq. 24 


Combining Eqs. 32 and 34, it follows that 


ro + qe < 2. 


From Egq. 35 several possibilities are suggested. If 7. = 1 and 


g: = 0, then Eqs. 32 and 33 become 


< 3. (36) 


If r. = O and q, = 1, then 


(3/) 


Neither Eq. 36 nor 37 can be satisfied by integral values of g; and 7, 
so that the assumed values of 72 and g2 are impossible. On the other 
hand, if r. = 0 and gz = 0, then 0 < g; — 7; < 1, which can be satis- 
fied if g. = 7, +1. Therefore, the final decomposition must have all 
double zeros = 0, r2 = 0) withg = r+ 2. Butif Ay(s) and By(s) 
have only double zeros, their sum can realize only complex zeros. 
All the real zeros of N(s), then, must have been deleted by Step 2 in 
the procedure. Therefore, the optimum decomposition has the form 


n m 


N(s) = II (Ss + (S Q 1)” + Ba, II bo, i)? ] (38) 


l 


where the a,-; and b,_, alternate since the optimum must satisfy 
Theorem I. Equation 38 is of the same form as Eqs. 20 and 21 so the 
theorem is proved. 


THE EXISTENCE OF THE OPTIMUM 


The existence of the optimum decomposition will now be established 
by showing that it may be simply formed from N(s). Since the real 
zeros of N(s) are separate from the remainder of the decomposition 
(Eq. 38), the only question is whether the complex zeros may be de- 
composed as shown. 


Theorem III: Any polynomial of the form 


m 


N(s) = JI (s + s,)(s + &,) (39) 


1 


J. 
2 + — — fe (33) 
> — 71) + G2 — fre. (34) 
(35) 
m—1 | 
l 1 l 
1 
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can be decomposed into the form 


N(s) = Ai(s) + Boi B,(s) (40) 
where 
Ax(s) = J] (s + ax-1)? (41) 
and 
m—1 
By, By(s) = Bo, II (s + 1)? (42) 
1 
where 
1. Bo, is positive and real, (43) 
2. dox—1 and are real, (44) 
3. < < a3 < 03 < +++ < (45) 
Proof: Let 
A,(s) = [Rn»(s) (46) 
and 
By:By(s) = [On(s) ? (47) 
so that 
R,,(s) = Il (s + (48) 
1 
and 
m—1 
Q,(s) = VBor I] (s + (49) 
1 
Then 
N(s) = [Ra(s) + [Qm(s) (50) 
= [Rn(s) + jQm(s) ][Rm(s) — jQOn(s) ]. (51) 


Now consider N(s). Its factors are of the form 
N(s) = (s + jd) (s +c: jdi)] (52) 
1 1 


where c; are real and d; > 0. Comparing Eqs. 51 and 52, the following 
identifications can be made 


R,.(s) + jOn(s) = (s tes + jd.) (53) 


é 
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m 


— 7Qm(s) = TI (s +e; — (54) 


1 


It must now be shown that the coefficients and zeros of R,,(s) and 
Q,,(s) satisfy Conditions 43-45. 
First consider the case m = 2. 


R2(s) + jQ2(s) = (s + jdi)(s + ¢2 + jd) 
= (s+ c¢,)(s + — did 
+ jldi(s + cs) + +c,)]. (56) 


Hence, 
R.(s) = (s + €1)(s + €2) — didz (57) 


and 


Q2(s) = Co) d»(s C1). 


root locus of R.(s) and Q2(s) as d; and d», are varied is shown in 
Fig. 5. 


Root Locus of R_(s) 


Root Locus of Q, (s) 


Fic. 5. 


The coefficients of the highest power of s in R2(s) and Q2(s) are 1 
and d,; + ds», respectively, both of which are positive and real. More- 
over, the zeros of R.(s) and Q2(s) are shown by Fig. 5 to be real and to 
alternate along the real axis. ‘Therefore the theorem is proved for 
m = 2. 

For m > 2, suppose, in the process of forming R,,(s) and Q,,(s) by 
multiplication of terms like s + c; + jd;, the r** multiplication is being 
performed. Then 


r—2 


R,(s) + jQ-(s) = (s ter + jd) (s +ai) + jK,-1 II (s + (59) 


1 


120 
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R,(s) + 70,(s) = [(s + TI (s + a,;) — d,K,_, Tl (s + B:) ] 


r—1 


+ fla, Th + TG + 


r r—1 
, 
= TJ (s + jK, TI (s + B/’). (61) 
1 1 
j¥ 
Cr 
Root Locus of 
Cr 
Root Locus of 
Root Locus of ~Bp $C, S “Op 
-Cr 
Root. Locus of Q (s) 
(b) 
-Cr 
4-0, << 
Root Locus of Op Cr Bp-1 
- “Cr Ss “Bp-1 


Root Locus of Q (s) 
r 
(c) 


Fic. 6. 
It will now be shown that, if K,_1, a:, and 8; satisfy conditions equivalent 


to Conditions 43-45 then K,, a;’, and 6;’ do also. 
Consider the root locus equations of R,(s) and Q,(s): 


r—2 
d,K, 1 II (s + B:) 
1 


(62) 
(s + ¢,) II (s + aj) (s + ¢,) II (s + aj) 
| +e) (s 
— =14— 


d, II (s + ai) d, II (s + a) 


Now d, > 0 and K,_; > 0 by assumption that the latter satisfied Condi- 
tion 43. Since the alternation of a; and 6; along the real axis is also 


: 
1 
r—2 
1 1 
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assumed, the form of the root loci of R,(s) and Q,(s) can be shown 
(Fig. 6). In Fig. 6a, the zero at —c,, which may lie anywhere on the 
real axis, is assumed at infinity. Clearly the zeros of R,(s) and Q,(s) 
alternate in this case. If now the zero at —c, is brought from infinity 
along the negative real axis, it will not affect the permitted regions of 
“closed loop” zeros to the right of it. To the left of it, —c, will reverse 
the permitted locations of the zeros of R,(s) and Q,(s). Therefore, 
alternation of the zeros of R,(s) and Q,(s) is preserved in the presence of 
—c, That the zeros of R,(s) and Q,(s) remain real is shown by Figs. 6) 
and 6c, when —c, is located to the right or left of any —a,. 

The theorem has now been proved first for ry = 2 and then for 
2<r<m. Therefore, it is certain that at each multiplication by a 
term like (s + c; + jd;), the result will satisfy conditions equivalent to 
Conditions 43-45. If multiplication is continued, then, until 7 = m, 
the final By;, —a;, and —b; fulfill Conditions 43-45. Therefore, the 
theorem is proved. 


THE UNIQUENESS OF THE OPTIMUM 


If more than one decomposition exists of the form of the optimum, 
then it must be possible to form one such decomposition from the other 
by the addition or subtraction of a polynomial. Then 


Avyii(S) = Ay(s) + 605(s) (64) 


Bons iBy, i(S) = BonBy(s) 60d (S). (65) 


Since A y(s), By(s), Awai(s), and By,,(s) have the real zeros of N(s) in 
common, 6(s) must have these zeros also. These zeros cancel in the 
above equations and may be ignored from further consideration. The 
equations may then be written 


m m 


II (Ss + a's 1)” II 1)” 600(S) (66) 
1 1 


ll 


m—1 m—1 
Bows II (§ 4- b’», 1)" Bown II 4- 1)" 600(S). (67) 
1 1 


The sign preceding 6) may be fixed since the decomposition identified 
with the left hand side of these equations is arbitrary and may be identi- 
fied as the one for which the above sign choice is correct. If 


5(s) = JI (s + 6,) (68) 
1 


then t < 2m — 2 or else Boy,, will be negative. Now if the root locus 
of Eq. 67 is examined (Fig. 7a) as 69 varies, the regions along the nega- 
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~ A s 
(a) Root Locus of and 


© 


(b) Regions Permitted for Zeros 


Fic. 7. 


tive axis permitted for the zeros at —b’.,_; may be determined solely 
from the zeros of 6(s). This is due to the duplicity of all the ‘“‘open- 
loop” zeros at —b,_;. For if the root locus condition 


are n= 0, +2, +4, (69) 


IL (s + 
1 


is satisfied to the immediate right of —b.._,; then it will be satisfied to 
the immediate left of —b,_;. Likewise, if Eq. 69 is not satisfied to the 
immediate right of —b»_;, then neither will it be satisfied to the left. 
Similarly, the regions permitted for —a’s,; are determined from the 
satisfaction of the equation 

arg = n= +3, +5,---. (70) 


m 


(Ss + 


Again, the satisfaction of this equation on the negative real axis depends 
solely on the zeros of 6(s). However, since the values of ” that may 
satisfy Eqs. 69 and 70 are different, the regions permitted for the zeros 
at —a’»,; and —b’.,_; do not overlap. 
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Now assume that all the zeros of 6(s) are real and distinct and that 
there are 2m — 2 of these zeros. In Fig. 7), the zeros permitted in 
each region are shown in brackets. Note that there is room for only 
2m — 2 zeros of Aw4i(s). If any zeros of 6(s) are multiple or complex, 
there are even fewer vacancies for the zeros of Ay,;(s). The conclu- 
sion is that no 6(s) exists which permits formation of an Ay,,(s) of the 
form shown. The uniqueness of the optimum decomposition follows. 


THE ANGLE CRITERION 


To test whether a polynomial N(s) could be formed as in Eq. 2, 
A,(s) and Bo:B,(s) may be formed and —a,; examined. The non- 
positiveness of —a, is then a necessary condition that N(s) be so formed. 

Even so, there is an easier method of testing N(s) which is related 
to the manner in which the optimum decomposition is formed. In 
addition, this different testing procedure is interesting in the simplicity 
and compactness of the demand it imposes on N(s). The equivalence 
of the two testing methods is shown by the following theorem. 


Theorem IV: For any N(s) of the form of Eq. 39, a necessary and 
sufficient condition that —a; be non-positive is‘ 


m 
LX arg si < 5. (71) 


1 


= c¢; + jd;, define 


Proof: For any s; 


wi(o) = Vd? + (ci + 0)? 


= tan-! — 0 < <7. 


Ci: +o 


The graphical interpretation of the above quantities is shown in Fig. 8a. 
Now define 


= (74) 


so that it must now be shown that 


args; = 0,(0) = y(0) < (75) 
1 


1 


4 This condition is identical with the results of Horowitz (3) for a fourth order polynomial. 


te 
(72) 
: 
: 
(73) 
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Let the s-plane be shifted by an amount ¢, and form the optimum de- 
composition in the shifted plane. 


Ii[s +(e) (77) 


Interpretation of (0), (0) 
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The constant term of R,,(s + ¢) is 


r(o) = [J (ax_-1 + ¢) = Il w,(o) Cos 


m 


= w,(c) cos y(o). 
Therefore 


r(a) 0 
(79) 


Now y(oe), from Fig. 8a, is a monotonically decreasing function of o, 
achieving its maximum of mz at o = — « and its minimum of zero 
at 0 = + ». Therefore, if y(o) is plotted (Fig. 85) beginning at 
o = + «, then clearly y(—a,;) = 7/2. Hence, if —a, < 0 


MI —(1,) 


On the other hand, if y(Q) < 7/2 
y(O) < 2/2 = y(—a,) (81) 


which implies —a, < 0 from the monotoneity of y(¢). Therefore, 


and the theorem is proved. 
REALIZATION OF ONLY COMPLEX NATURAL FREQUENCIES 
When a polynomial N(s) has only complex zeros and these satisfy 
the angle criterion, the optimum decomposition may be applied to 
yield 


m 


m—1 
N(s) = II (s + an_1)? + Bo (s + (83) 
1 


The zeros of Ai(s) and B,(s) are readily manipulated to form 


m—1 
VBor IL (s + 
Zrec(s) = — — 


II (S + 


= 

Be 

2 
: 
Me 
(84) 
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and 


Il (S + Gex-1) 


TI] (s + 
1 


which are RC and RL driving point impedance functions. 


REALIZATION OF REAL AND COMPLEX NATURAL FREQUENCIES 


The optimum decompositon of an N(s), the complex zeros of which 
satisfy the angle criterion, is, according to Theorem II, 


N(s)=II (s+p.) (s+p,) TL (86) 
1 1 1 1 


The problem of manipulating the real zeros of N(s) now becomes 
apparent. They cannot be assigned to the same impedance, as for 
example 


n m—1 
Il (Ss -+}- Pi) Il (s be, 1) 
1 


m 


n (87) 
II (s + pi) Il (S + dox—1) 


II (s + 


II (s + 
1 


for they obviously cancel in Eq. 87. If, instead, these zeros are as- 
signed to each impedance 


n m—1 
I (s + pi) I (s + bex-1) 


Il (S + Ax-1) 


II (s + pi) Il (s + ax 1) 
II (s + 


then Zrc(s) and zrz(s) cannot both be RC and RL impedance functions, 
respectively. 

Since it seems likely that there is some other decompositon of N(s) 
that satisfies Theorem I and that also can be manipulated to form a 
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Zrc(S) and a Zprz(s), the question arises as to how this decomposition 
can be formed from the optimum. There are many ways of forming such 
a decomposition by choosing various 695(s) in the equations 


A onAn(s) = Ai(s) — 65(s) (91) 
= BoBi(s) + (92) 


One such method of selecting 64(s) is the following: 


1. Locate a zero of 6(s) at every point where 


S = — and 

2. Locate a zero of 6(s) in the immediate vicinity of every s = — pi. 
These zeros are chosen to the right of, left of, or else at s = — p; in 
such a manner that the resulting Ay(s) and By(s) satisfy Theorem I. 
Figure 9a shows an example of incorrect choice of the zero at s = — 5s. 
Theorem I is then violated at s = — o,; for the resulting Ay(s) and 


By(s). A correct choice of the zeros of 6(s) is shown in Fig. 9b. Should 
there be a zero of order / at —p; (— pz» in Fig. 9c), locate / — 1 zeros of 
5(s) at this point and proceed to decompose the remaining zeros of 
A,(s) and B,(s) according to the above steps (Figs. 9d and 9e). These 
/ — 1 zeros then are deleted from the root locus of Eqs. 91 and 92 since 
they are contained in A,(s) and B,(s). After selection of a sufficiently 
small 69, By(s) and Ay(s) will be found to have / — 1 zeros at —p». 
To separate these / — 1 zeros, form a new decomposition from A y(s) 
and By(s) just as the latter were formed from Ai(s) and B,(s), respec- 
tively. The 6(s) chosen is the same as the original except that a zero 
of 6(s) is deleted at —p. (Figs. 9f and 9g). The root loci of this step 
in the region of — 2 are shown in Fig. 94. It should be noted that the 
zeros of Ay(s) and By(s) uncovered by the deleted zero are seen in 
Fig. 9h to move in the same direction. However, the form of the re- 
sultant decomposition must be that of Fig. 97 to realize the double 
zero of N(s) at —p2. The decomposition of Fig. 97 cannot realize a 
double zero at that point. This step of uncovering a zero at s = — p» 
may be repeated again and the zeros of the new decomposition in the 
region of —.2 are shown in Fig. 9k. Repeating until the /* zero is 
separated, the distribution becomes that of Fig. 9m if / is odd and that 
of Fig. 9n if liseven. In the former case, there are / — 1 pairs of zeros 
whereas the latter results in / pairs of zeros in the region of —p2. The 
philosophy in performing these repeated decomposition formations is 
that a sufficiently small 69 is chosen each time such that no zeros of 
Ay(s) or By(s) become double and such that no zeros of Ay(s) enter 
the right half plane. This is clearly always possible since any 659 > 0 
will achieve separation of the real zero at s = — py. 
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The resultant Ay(s) must always, regardless of 65(s), have a zero 

more positive than the most positive zero of Ai(s). Theorem II has 

shown that Ay(s) can always be reduced to A,(s) while decreasing the 

most positive zero of Ay(s). Therefore, if Ai(s) has a zero at s = 0, 

A y(s) must contain a zero in the right half plane. It is necessary, then, 
m 


that > arg s; < 2/2 and also that all the real zeros of N(s) be negative 


1 
in order to achieve an impedance decomposition. These conditions 
are sufficient as well as necessary since the optimum decomposition of 
any N(s) satisfying them can be converted to an acceptable decomposi- 
tion by choice of an arbitrarily small 6) in Eqs. 91 and 92. 


An Example 


The polynomial 
N(s) = (s? + 2s + 2)(s? + v3s + 1) (93) 


is to be formed as in Eq. 1. 
Separating the upper half plane zeros, one identifies 


R(s) + jQ(s) (s + 1+ j1)(s + 867 + 7.5) (94) 
(s + 5:)(S + 52). (95) 


Here, arg s; + arg s2 = 52/12 so that RC-RL decomposition is pos- 
sible. Further, 


R(s) + jQ(s) = (s + .226)(s + 1.64) + j1.5(s + .91) 


(96) 


so that 
A,(s) + BoiBi(s) = (s + .226)?(s + 1.64)? + 2.25(s + .91)?. 
The desired formation of N(s) is then 

N(s) _ (s + .226)(s + 1.64) 1.5(s + .91) 


D(s) 


(97) 


1.5(s + .91) 


(s + .226)(s + 1.64) 


CONCLUSIONS 


In summary, the necessary and sufficient conditions that a given 
polynomial V(s) be formed from the sum of RC and RL driving point 
impedances are stated: 


(99) 


1. If N(s) = (s s,)(s + then ate si < 
1 1 


2. If N(s) Il (s + pi) Il (s + s,)(s + 8,), then: 
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(100) 


a. JJ args; <3 
1 


b. pi > 0. (101) 


The sufficiency of the above conditions may be expected to have appli- 
cation in both passive and active network synthesis, where the parti- 
tioning of network functions frequently leads to simplified models. 
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SPACE, TIME, AND ETERNITY 
BY 


GUSTAF STROMBERG '! 


According to the theory of relativity, measurements of space and 
time intervals (determined with yard sticks and clocks) are dependent 
on one another in such a way that, mathematically, they can be repre- 
sented as being ‘‘projections’’ of a four-dimensional world-frame called 
space-time. This term was introduced in 1908 by the famous German 
physicist H. Minkowski, who pictured space-time as a ‘“‘ghostly”’ 
reference frame in an external world in which all physical phenomena 
were displayed. According to Kant, space and time are a priori forms 
of human perceptions and are apparently the product of the sensory 
nerve system in the human body. Of the greatest importance for our 
knowledge of the world around us are our sense of vision and our sense 
of the flow and duration of time which have enabled us to form ideas 
about the existence of external objects and of their motions. The 
science of physics has been developed for the purpose of co-ordinating 
our space and time measurements so as to furnish a logically consistent 
set of rules, generally called ‘“‘laws of nature,’ governing these measure- 
ments. Minkowski’s four-dimensional space-time was used with great 
success by Einstein in his general theory of relativity, to formulate laws 
of motion of greater generality and exactness than those formulated by 
Newton. 

But, there are many phenomena in our experience which cannot be 
expressed in terms of measurements of space and time intervals. Be- 
cause they are not governed by the laws of physics, such phenomena may 
be called nonphysical. A simple example is the moving pictures we 
see in a dream or an hallucination. Because of their private and per- 
sonal nature they are classified as subjective, and are therefore regarded 
as having no objective reality. Such pictures, however, are built of 
the same kind of ‘‘mental stuff” as are the pictures we see with our eyes 
when we are fully awake and not under the influence of hypnotic sug- 
gestion or drugs. 

In spite of the great advance in our knowledge of the physical 
world of space and time in recent years, little progress has been made in 
the field of research which deals with the relationship between the physi- 
cal and the nonphysical worlds, sometimes referred to as the relation- 
ship between mind and matter, or, in popular language, between our 
body and our “‘soul.”’ 


! Retired member of the scientific research staff of the Mt. Wilson Observatory of the 
Carnegie Institution of Washington; now in Pasadena, Calif. 
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The world of matter and the world of mind appear at first sight so 
different from one another that it seems impossible to find any con- 
necting links between them. In attempting to explain ‘heir relation- 
ship, some philosophers have tried to describe consciousness as an 
“epiphenomenon,” that is, as a by-product associated with the stimula- 
tion of sensitive nerve cells in the cortex of the brain. Some philoso- 
phers claim that what we call ‘“‘the physical world” is a play of mental 
““shadows,”’ seen directly or pictured in our imagination, which we must 
be careful not to identify with ultimate reality, which itself may well 
be completely beyond human comprehension. Some scientists claim 
that the rules governing the mental shadow-play represent a more 
fundamental level of truth than do the individual sense data. 

To represent the relationship between quantities of different nature, 
it is often convenient to plot them as co-ordinates in a reference frame 
with axes perpendicular to one another. For instance, the path of a 
moving object in the space-time world can be pictured as a continuous 
line in a four-dimensional reference frame. The point where two such 
“world lines’’ intersect represents an observable ‘‘event”’ in the physical 
world of space and time. A mental event, on the other hand, takes 
place at a certain time, but not at any particular location. Therefore, 
to include mental events in our world picture, we must add to the four 
dimensions of space-time a new “‘dimension”’ or ‘‘axis’’ perpendicular 
to the three space axes and to the time axis. (Mathematically, there 
is no difficulty in having as many independent dimensions in our world 
picture as we think are needed.) 

The idea of a five-dimensional universe has been found very useful 
in theoretical physics. The German physicist A. Kaluza (1)? showed 
in 1921 that a unified field theory could be developed if a fifth dimension 
were added to the four-dimensional world of Minkowski and Einstein, 
and similar conclusions have been reached by other theoretical physi- 
cists. Of special interest is the fact that Prince Louis de Broglie, the 
prominent French physicist who in 1924 introduced the idea that all 
material particles must be associated with some kind of “waves,” 
called attention to the fact that the hypothesis of a five-dimensional 
world-frame was a logical generalization of the relativity theory, and 
that it involved the assumption that the human sense organs are in- 
capable of directly perceiving any effect due to the existence of the fifth 
dimension, and that what we actually perceive must be the “‘projec- 
tions”’ onto the space-time world of events occurring in an extra-physical 
world, that is, a world beyond physical space and physical time (2). 

An important step in the change-over from a four-dimensional to a 
five-dimensional universe was taken by the prominent British scientist 
and philosopher, J. G. Bennett. In collaboration with two British 
mathematicians, R. L. Brown and M. W. Thring, he published in the 


2 The boldface numbers in parentheses refer to the references appended to this paper. 
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Proceedings of the Royal Society for July, 1949, an article in which he 
showed that the mathematical foundations of theoretical physics could 
be simplified and generalized by the assumption that the universe had 
five dimensions.* 

Bennett calls the new axis the efernity axis, and it runs perpendicular 
to Minkowski’s time axis; since it is also perpendicular to the three 
space axes, it defines a realm beyond the four-dimensional space-time 
world of physics. We may call this recently discovered world the 
“Eternity Domain of the Universe.”’ Bennett tells us that the reason 
why we have no immediate experience of events which take place in this 
nonphysical realm is that we are “eternity blind.” 

Certain properties of the eternity domain can be deduced from the 
fact that it lies beyond both space and time in their physical sense. 
No physical measurements can be made in that world, because measure- 
ments are made with material measuring rods and clocks, and matter 
has no place in a nonphysical world. Since there can be no units of 
length or of time in the eternity domain, it must be described as non- 
metrical. Separations in space cannot be defined, although there may 
be some kind of ‘“‘mental separation.’’ Long and short time intervals 
cannot be distinguished from one another, so that duration in time is 
meaningless. On the other hand, there are reasons to believe that a 
sequence in time has a well-defined meaning (see below). We may expect 
that there is a causal relationship between events in the physical and 
the nonphysical worlds. 

We shall now make a fundamental assumption. In the nonphysical 
world beyond space and time, which has here been called the eternity 
domain, lies the ultimate origin of all things: energy, matter, life, 
consciousness, and mind. In short, all the characteristics of the world, 
physical as well as mental, are supposed to have their origin and their 
sustaining “‘roots’’ in this recently discovered extra-physical domain. 
This assumption should be regarded as a “working hypothesis,’’ and its 
justification depends on whether or not it can explain, or rather help 
us to visualize, the relationship between physical and mental phenomena. 
In the following pages are given some examples of how the introduction 
of the eternity domain can help us to picture the nature of this 
relationship. 


1. THE LOCATION OF MEMORIES 


The following questions have often been asked: How is it possible 
for us to remember events which happened long ago? Where and in 
what form are our memories stored, and how are they brought forth 


3 The story of the development of the theory of a five-dimensional universe, and how 
Bennett was influenced by the ideas of the Georgian mystic, Gurdjieff, is related in Bennett's 
book, “The Crisis in Human Affairs,’’ London, Hodder and Stoughton, 1948. An extensive and 
technical description of the new ideas is presented in Bennett’s more recent book, ‘‘The Dra- 
matic Universe,”” London, Hodder and Stoughton, 1956. See also Maurice Nicoll’s ‘Living 
Time,” published by Vincent Stuart, London, 1952. 
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in our conscious minds? It has sometimes been suggested that our 
memories are recorded in some form in the material brain, and that they 
are evoked by the association of ideas, by other memories, by visual 
observations, or, as recently found, by direct electrical stimulation of 
certain nerve cells in the cortex of the brain. But the. brain is a fluid 
mass and must therefore, except at extremely low temperatures, be 
subject to the disruptive effects of the thermal motions of its molecules. 
Further, it is impossible for us to understand how the material structure 
in the brain can produce mental effects in our consciousness. 

The introduction of the eternity domain has given us a clue to the 
solution of this difficult problem. We may think of the brain, as re- 
vealed by our sense organs and by our scientific instruments, as being a 
“projection” onto the space-time world of its inherently nonphysical 
counterpart in the eternity domain. We must always remember that, 
in whatever way our knowledge of the physical brain has been acquired 
this knowledge is in itself an zdea of something which we can neither 
understand nor adequately describe. What we see with our eyes can 
only be a partial aspect of an underlying reality. 

If a man’s memory complex belongs to the immaterial eternity do- 
main it can be expected to be free of any deterioration. This can ex- 
plain the well-known fact that memories from our childhood can persist 
during a long life, in spite of the elimination of waste products, which 
goes on continuously during the metabolism in all living structures, 
so that we have a completely “new” brain after a few years. 

The evocation of a particular memory element can now be pictured 
in the following way. A group of nerve cells in the brain is ‘‘stimulated”’ 
or ‘‘activated”” by some means, for instance by being supplied with an 
extra amount of oxygen-carrying blood. The energy released in the 
oxidation or burning of the nerve tissue and the related modification 
of its electronic structure are associated with a corresponding change 
in the “spiritual’”” memory complex of the individual in the eternity 
domain (3, 4). The potential connections between the particular 
parts of the material brain and their spiritual counterparts then become 
effective, and the memory of a particular event appears in the conscious- 
ness of the owner of the brain. This idea is similar to Henri Bergson’s 
theory of memories (5)—that our memories ‘“‘are hidden from our 
knowledge by a veil of atoms.’’ Because the physical brain is built of 
material elements having inertia, it takes energy and time to stimulate 
a memory element, but when it is evoked, it persists for a certain time. 
The process of remembering and of thinking is more or less “‘sluggish,”’ 
so far as the brain mechanism is concerned, but in the nonphysical 
world these processes can be expected to be instantaneous. 

The theory here expressed is supported by evidence offered by the so- 
called panoramic death vision (6). It has been reported for many cen- 
turies that a person near death, for instance by drowning, sometimes 
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can see a panoramic view of his life, in some cases from childhood on. 
The moving pictures often show details which the person having the 
vision has forgotten long ago. It seems that when the oxygen supply 
to the brain has been cut off, and when the immaterial brain field is in 
the process of disentangling itself from its load of incorporated inert 
matter, normal mental sluggishness disappears. It is important to note 
that to an outsider the vision seems to last only a very short time, while 
to the participant, who later tells the story of his vision, it may have ap- 
peared to represent a long time interval, when a great number of 
things happened. Such death visions may be much more frequent 
than we think, since reports of such visions are obtained only from those 
who are resuscitated, and not from those who perish immediately 
following such visions. It is highly significant that the pictures seen 
are quite clear and do not overlap, as in a well-tuned television instru- 
ment. It indicates that in the eternity domain a sequence in time is 
meaningful, but a duration is not. 


2. THE ORIGIN OF MATTER AND ENERGY 


It is generally admitted that present-day science cannot give us any 
information about the ultimate origin of the matter and energy in the 
world of today. For instance, Hoyle’s theory of a continuous creation 
(7) and the theory of a cosmic explosion developed by Alpher and Her- 
man (8) give us no hint about the ultimate origin of energy. The intro- 
duction of an eternity domain in the universe makes it possible for us 
to form a mental picture of the origin of energy, and consequently of 
matter. 

Let us think of the eternity domain as a continuum everywhere in 
contact with the four-dimensional space-time world, somewhat like the 
atmosphere in contact with the surface of the earth. At certain contact 
points, which may be called energy sources (9), energy emerges from the 
eternity domain into the space-time world. This energy may be thought 
of as previously existing in an undifferentiated and potential form, which 
can best be described as mental. At the very moment of their emer- 
gence, the elements of energy probably appeared in an extremely con- 
centrated form, possibly as neutrons, as in Alpher’s theory. Before 
their emergence, their location and size were undefinable, because these 
concepts have no meaning in the eternity domain, which lies beyond 
physical space. In the theory of continuous creation, the emergence 
goes on at all times, while in the explosion theory it took place at a 
particular time and with extreme violence. This explosion, if it actually 
occurred, represented the “birth” of the physical world and the begin- 
ning of physical or clock time. Recent studies indicate that this be- 
ginning of time took place between four and five billions of years 
ago (10). 
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3. THE ORIGIN OF LIFE AND OF CONSCIOUSNESS 


Most biologists maintain that organic life on the earth is the result 
of a slow development which started with the formation of certain 
chemical compounds in the fluid water at the surface of the earth, the 
needed energy presumably being supplied by sunlight. From these 
simple beginnings, living cells were formed in some unexplained way 
and from them came higher forms of life, such as plants and animals. 
In the struggle for survival, a natural selection of living forms took 
place which favored the strong and hardy individuals, while the weaker 
and unfit forms were eliminated. This development was made possible 
by the close similarity between parents and offspring, coupled with 
seemingly accidental changes in their gene systems, the carriers of the 
hereditary factors. 

In 1935 Burr and Northrop of Yale University proposed the “‘electro- 
dynamic theory of life’’ (11). This led to an extensive investigation by 
Burr and his associates at the Section of Neuro-Anatomy of the Medical 
School of Yale University of the fine structure of the electric fields in 
which all living plants and animals are known to be imbedded. Many 
important facts were discovered during this investigation (12, 13). It 
was found that the electric pattern had a complex structure which 
extended well beyond the visual limits of the living organism. If the 
oxygen supply necessary to maintain the metabolism was reduced, the 
field contracted, without change in structure, and disappeared com- 
pletely at the death of the organism. Burr stated (13) that “‘it is hard 
to escape the conclusion that the electrical pattern is primary, and, in 
some measure at least, determines the morphological pattern.”’ It 
seems that during growth fluid matter in the form of nourishment of 
particular electric structure (for example, proteins) is incorporated in 
a pre-existing electric field of complex structure and that it becomes 
directly observable by this form of “materialization.” 

Observations of this sort led the author in 1945 to propose the theory 
of ‘‘autonomous, living fields’ (14), based on the assumption that the 
sources of the organizing electric fields cannot be located in the 
electrically charged particles of which matter is generally supposed to 
be built up, but should be regarded as singularities in a pre-existing 
force field with topological or morphological, but with no metrical, 
properties, such as size. Fields of this sort can be imagined as existing 
in an extremely contracted and dormant form in an egg cell or seed. 
Going one step further it may be postulated that the living field, when 
it is in its potential form, has no size at all, in which case it must be 
regarded as a nonphysical potentiality.‘ 

The idea of the existence of ‘living sources’ as the points of entrance 


4 Bennett, who regards this as a sixth dimension in the universe, calls it Hyparxis, which 


means “‘ableness to become.” 
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in full harmony with the theory of ‘emergent energy” (9). Life seems 
to have emerged from ‘‘another world” than that described in the science 
of physics. The idea is almost identical with that presented by 
Northrop and Burr, in which the organizing fields are ‘‘primary proper- 
ties of the universe, and that they are, to a greater or less extent, 
modified by the presence of matter, the field and the particles being 
thus dependently related.’’ (Abstract of a letter from Burr to the 
author.) 

The assumption that life originates in the eternity domain can help 
us to understand two important problems: the seemingly purposive 
activities in living organisms, and the ultimate origin of consciousness 
and mind. These two problems are intimately related, I believe. 

The theory that the growth and development of a living organism is 
directed by a purpose is called teleology, which can be translated as 
“looking ahead”’; it is particularly evident during embryonic develop- 
ment. Teleology can be regarded as the opposite to mechanism. With 
the general acceptance of mechanistic ideas, teleology has generally 
been regarded with much skepticism and often with hostility by most 
biologists, since it may open the door to mysticism of one sort or another. 
But some prominent biologists have accepted some form of teleology, 
for example, E. G. Conklin (15), who says ‘‘There may be good grounds 
for holding that mechanism and teleology are complementary views of 
nature, neither excluding the other.’” The problem of randomness 
versus directiveness has been discussed in great detail by Ralph S. 
Lillie (16) and similar ideas have been presented by E. W. Sinnott (17). 

There is no reason to believe that an aggregation of atoms and 
molecules, simply because of its chemical structure, however complex 
this structure may be, can acquire any form of consciousness, not even 
of the simplest form. Asan example, take the ability of men and many 
animals to perceive light and colors. It is known that this ability is a 
consequence of the existence of sensitive nerve cells in the retina of the 
eyes. From these nerve cells run nerve fibers to special neurones in the 
brain. When the nerve cells are stimulated by one means or another, 
be it by electromagnetic radiation, popularly known as “‘light,’”” by a 
hard blow on the eyeball, by chemical irritants, or by static electricity 
applied directly to the optic nerve cells in the brain, there results a 
sensation of light or colors. How can the causal relationship between 
the mechanical stimulation and the subsequent mental sensation of 
light be explained? The hitherto unsolved problem of the relationship 
between mind and matter appears here in all its poignancy. 

The introduction of the eternity domain, in which consciousness has 
its origin and its sustaining ‘‘roots,”’ may help clarify the nature of this 
relationship. 

A particular type of “living elements,’’ which may be termed 
“potential optic nerve cells,’’ have emerged from the mental eternity 
domain into the physical space-time world. After they became associ- 
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ated with an immaterial nerve system, and after the corresponding living 
fields had expanded and become materialized by the incorporation of 
matter of proper structure, they became parts of the gene system in the 
egg cell of the animal. (The fact that color blindness is linked to sex, 
indicates that certain “color genes” are located in the X-chromosomes. ) 
When fully developed, the optic nerve cells become located in the retina 
of the eyes of the animal and are connected with complementary nerve 
cells in the occipital part of the brain. The connection between the 
optic nerve cells and their “cosmic” or “‘spiritual’’ counterparts in the 
eternity domain had never actually been lost, and when the nerve cells 
were stimulated by one means or another, which involves an absorption 
of energy, a ‘“‘contact” between the physical nerve cell and its spiritual 
origin was established. As long as the field was in its activated and 
therefore sensitive state, a pre-existing potential association between 
the physical and the nonphysical world had become an actuality, and the 
animal or the man perceived a sensation of light or colors. Since there 
is reason to believe that the eternity domain forms an integrated and 
nondivisible whole, it may be said that a particular color sensation has 
its origin and its ‘‘roots” in the World Soul,’ which itself is spiritual and 
purely mental. (This explanation of the origin of color sensations was 
first presented in the author’s book, ‘“The Soul of the Universe” (18).) 


4. THE IMMORTALITY OF THE HUMAN SOUL 


Peoples of many races and religions have expressed belief in the 
immortality of the human soul. The origin of this idea may lie in the 
hope of again seeing departed friends and relatives; it may have been 
fostered by assertions by wise or “holy’’ men; or it might have been due 
to a subconscious knowledge of the existence of another world than 
that experienced by our sense organs. It has been reported from many 
sources that messages have been received from the “spirit world”’ which 
contain information known to no one but the departed person, the in- 
formation being later verified by an investigation. Natural science has 
reacted unfavorably to such ideas, largely because present-day science 
can give no explanation of such phenomena. 

The introduction of an extra-physical dimension in our world picture 
has made it possible to form a scientific picture of the survival of the 
human soul at death. When the immaterial force field in the brain 
contracts, as do other living fields in the absence of oxygen, the contrac- 
tion proceeds so that this field finally has no size at all, although its 
essential properties remain intact. When it has completely disappeared 
from the physical world, it becomes what may be described as a “‘living 
source.”’ In other words, zt returns to the mental world from which tt 
originally came. As evidenced from the ‘“‘death visions’ described 
earlier, all memories in their minutest details can be recalled in an 


5 The idea of a World Soul will be explained later. 
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instant. This indestructible memory complex forms an essential part 
of the indivisible entity we call a soul, which apparently does not have 
to be associated with matter. It seems that the eternity domain is 
our real “home,” and from there we make short excursions to the physi- 
cal world of space and time. During such a trip or trips we learn many 
lessons which may be helpful in our further development. 


5. THE MECHANISM OF MENTAL TELEPATHY 


Thanks to systematic investigations, mainly carried on by Societies 
for Psychical Research in several countries and by the Psychology 
Departments in some Universities, particularly at Duke University, 
some facts have been established which make it necessary to postulate 
the existence of an extra-physical dimension, or rather a domain where 
physical separation in space seemingly has no effect. Experiments 
have shown that ‘“‘thought pictures’”’ can be transmitted over extremely 
long distances, and that the intensity and clearness of the mental im- 
pressions do not depend upon the distance between the sender and the 
recipient of the mental message (19). If a message were broadcast in 
all directions, and if it, while in transit, could be regarded as a particular 
form of energy, which could be described as “thought waves,” its in- 
tensity would be expected to decrease with the square of the distance 
from the transmitting person. Also, due to the interference of other 
thought waves, the reception should be blurred. The specific fact, 
verified by many hundreds of testified observations, that the image 
of a dying person can appear clearly in the mind of another person 
indicates that, under unusual conditions, a mental contact between two 
persons widely separated in space can be established. No evidence of 
any time-lag has been established. Such conditions lead to the con- 
clusion that the message or the mental picture does not “‘travel’’ in 
physical space as a form of energy, but in the extra-physical domain, 
where physical separation has no meaning.* Consequently, it is logical 
to believe that thoughts and mental pictures are transmitted as such, 
that is, in their original mental form. This might have been expected, 
since the eternity domain is itself a mental world. 


6. THE ETERNITY DOMAIN AS A LIVING BEING 


It is obvious that the mental characteristics attributed to the eternity 
domain are those manifested in the human mind. The human mind 
is not simply a collection of mental elements, like sensations, feelings, 
thoughts, and memories, but it is a highly integrated unit, to which may 
be given the name Person, a word which literally means the sound 
emanating from a masked individual on a theatrical stage. A person is 
generally regarded as a being with some degree of intelligence, and the 


6 The problem of long-distance telepathy has been thoroughly investigated by Sir Hubert 
Wilkins and Harold Sherman (19). 
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“mask”? mentioned above may be considered something which hides 
him from the eyes of the public. This integrated mental unit called a 
person has often been called a soul, or its equivalent in the various 
languages spoken on the earth. 

It has long been customary to speak about man as a combination 
of a soul and a body, both of which are highly integrated units. The 
integrated character of the human body is evidenced by the co-ordinated 
activities in our bodies, activities which reveal man not as a colony or 
aggregation of cells, but as a highly organized system in which the vari- 
ous parts perform specific functions necessary for the life of the organism 
as a whole. 

The eternity domain is regarded as the realm in which all mental 
characteristics are ‘‘rooted.”’ It is very difficult to picture this domain 
as a loose combination of mental elements, like various types of sensa- 
tions, feelings of love and hate, pleasure and pain, of satisfaction and 
remorse. <A characteristic property of the eternity domain is the ab- 
sence of space separation, and therefore it should not be thought of as 
having any separated “parts.’’ A separation in time may exist which 
separates the past from the future, although there is some evidence that 
even this kind of separation may be bridged under conditions which 
would make pre-cognition possible. The many striking phenomena of 
teleological guidance in the living world, as opposed to the unguided 
phenomena in the nonliving world, may be thought of as evidence of a 
foresight and wisdom far surpassing our own, but reflected in a small 
measure in the intelligence and purposive activities of the human race. 
Among men, the purpose of an action is often consciously known, while 
among animals purposive activities are often governed by instincts, 
which represent behavior patterns acquired through the experiences of 
a great number of generations and presumably preserved in some form 
and duplicated in the germ plasm of the race. 

In view of what has here been said about the meaning of the term 
person, it is justifiable to apply this name to the nonphysical and mental 
realm called the eternity domain. It can now be described as an 
Almighty, Wise and Living Person, the Creator of all things, physical, 
mental, and spiritual. In our mind we have an “image” of this Person, 
and an idea of His existence and nature. We are created in His image, 
but although this image is incomplete and our idea of it distorted and 
colored by our own animal instincts, it is still an image of something 
greater than ourselves and of His whole creation. We can also under- 
stand that this Person in His wisdom may select one or more souls to 
carry important messages and admonitions to other human souls, while 
they are still associated with animal bodies and therefore have the ani- 
mal instincts of selfishness and lust for unrestricted power. Some of 
these messages we can also hear directly when we listen to the voice of 
the “Cosmic Conscience.”’ They tell us unequivocaliy that the essence 
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of Divine law can be expressed in the simple admonition : ‘Love ye one 
another!” 

For many centuries men of nearly all races have conceived of a 
Being or Spirit greater than themselves from which all things have 
originally come (20). It has been given various names, such as the 
Great Spirit, Jehovah, and Allah. Plato called it the World Soul 
(Anima Mundz1), and Giordano Bruno gave his life in defending Plato’s 
idea. In Vedantic idealism it has been called Parama-Brahman, a 
Sanskrit word meaning ‘““The ultimate spiritual essence”’ (21). 

In the English language it is called God. 
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NATIONAL BUREAU OF STANDARDS NEWS 


A STANDARD FOR MEASURING THE pH OF BLOOD 


To increase the accuracy of pH measurements in the physiologically 
important pH range 7 to 8, the National Bureau of Standards is recom- 
mending a pH standard! for determining the pH of blood. A phosphate 
buffer solution, the proposed standard has a pH of 7.414 at 25 C. and 
7.382 at 38 C. These standard pH values, accurate to within +0.002 
unit, will be useful to the research medical scientist, to the clinician, and 
to the physical scientist investigating acid-base behavior in the middle 
of the pH range. 

To discover the relationships that exist among physiological function, 
pathological condition, and pH, medical and biological laboratories 
have long studied the acid-base relationships in blood and other physio- 
logical fluids. Because many of these fluids are well buffered, the 
detection of any systematic variation of acid-base balance with physio- 
logical condition would be expected to require precise pH measuring 
equipment. Such equipment is commercially available in the form of 
the newer highly-sensitive pH meters with glass electrodes. However, 
these commercial instruments measure only a difference in pH and will 
yield useful pH data only if a suitable reference standard is available. 
To meet this need, V. E. Bower, M. Paabo, and R. G. Bates of the Bu- 
reau’s solution-chemistry laboratory developed the proposed standard. 
Their work results from a continuing program? to extend defined pH 
values through the pH scale and over a wide temperature range. 

The new standard reference solution is prepared by dissolving 1.179 
g (air weight) of potassium dihydrogen phosphate and 4.303 g (air 
weight) of disodium hydrogen phosphate in distilled water and diluting 
the solution to 1 liter at 25 C. The molalities of the two salts are, 
respectively, 0.008695 and 0.03043, and the ionic strength is 0.1. 

Standard pH values—which were derived from emf measurements— 
have been assigned to the reference phosphate solution at intervals of 
5 degrees from 0 to 50 C. To avoid the theoretical and practical diff- 
culties inherent in the estimation of liquid junction potentials, emf 
measurements of cells without liquid junction were used. 

An arbitrary, nonthermodynamic convention concerning single-ion 
activity must be adopted when defining pH. Because no one definition 


1 For further technical information, see ‘‘A Standard for the Measurement of the pH of 
Blood, and Other Physiological Media,”’ V. E. Bower, M. PAABOo, AND R. G. BATEs, J. Research 
NBS, Vol. 65A (Phys. and Chem.), p. 267 (May-June 1961). 

2 “Standardization of the pH Scale,”’ NBS Tech. News. Bull., Vol. 31, p. 138 (1947). 


145 


x 
| 
: 
: 
: 
: 
: 
ate 


146 NATIONAL BUREAU OF STANDARDS 


has been widely accepted, useful comparisons of highly precise pH data 
from different sources have been impossible until now. However, 
Bates and Guggenheim have suggested a simple, useful way’ of defining 
the activity coefficient of the chloride ion, which is needed to assign 
standard pH values. This definition has recently been endorsed by 
standardizing groups in the United States and Great Britian. Accord- 
ing to this proposal, the single-ion activity coefficient (yo) of the chlo- 
ride ion in dilute solution is so defined that it closely approximates the 
mean ionic activity coefficient of simple uni-univalent strong electrolytes 
in aqueous solution up to ionic strength (J) of about 0.1. The exact 
relationship adopted is: 


ANI 
= 
44-1,5v7 


where A is a constant of the Debye-Hiickel theory, the value of which 
depends on the temperature. 

With this new convention, standard pH values accurate to within 
+0.001 unit can be assigned to selected reference solutions with pH 
between 2 and 12. The ionic strength of these solutions should not 
exceed 0.1. 

The standard pH scale established by the Bureau several years ago 
is defined in terms of four fixed points which were obtained with standard 
reference solutions. These primary standards‘ were found to have a 
satisfactory degree of internal consistency.’ The values assigned to 
the new pH standard for blood are also consistent with these four 
standards. Together, these standards define a pH scale which, though 
conventional, is practical, convenient, and useful. 


FAST ELECTRON DOSE FROM A NARROW BEAM SOURCE 


The National Bureau of Standards has made a calculation of the 
penetration in air of fast electrons from a beam source. The work® 
was sponsored by the Atomic Energy Commission and was part of a 
program to develop methods for calculating the dose from a source 
geometry involving two spatial variables. Possible applications of such 
calculations include the determination of the spatial distribution of the 
dose in therapy involving electron beams and the estimation of the re- 
solving power of large single crystals in high-energy radiography. 

5 R. G. Bates AND E. A. GuGGeENuEIM, Pure and Applied Chemistry, Vol. 1, p. 163 (1960). 

‘The materials from which the standard solutions are prepared are described in NBS 
Circular 552, Standard Materials. ‘This Circular may be ordered for 35 cents from the Super- 
intendent of Documents, U. S. Government Printing Office, Washington 25, D. C. 

5R.G. Bates, Chimia, Vol. 14, p. 111 (1960). 

6 For further information, see ‘Calculated Energy Dissipation Distribution in Air by 
Fast Electrons from a Gun Source,” by Joun E. Crew, J. Research NBS, Vol. 65A (Phys. & 
Chem.), No. 2, p. 113 (1961). 
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The experimental situation is that of a collimated beam of mono- 
energetic electrons originating at a point in the medium. As the elec- 
trons move about, their energy is dissipated. The basic quantity to be 
determined is the energy dissipated in a small volume located at speci- 
fied distances from the source point and the center line of the initial 
beam. 

Calculations of the energy dissipation distributions for fast electrons 
in infinite and homogeneous media have been reported for the simplest 
source geometries.7"* The present calculation represents an extension 
of the earlier work and is a special application of the moment-fitting 
techniques developed at the Bureau and applied initially to problems of 
X-ray penetration.?:” 


7L. V. SpeNcER, Phys. Rev., Vol. 98, p. 1597 (1955). 

§8L. V. Spencer, NBS Mono. 1 (1959). For sale by the Superintendent of Documents, 
U. S. Government Printing Office, Washington 25, D. C. Price 45 cents. 

9L. V. SPENCER AND U. Fano, J. Research NBS, Vol. 46, p. 446 (1951). 

10 FANO, SPENCER, AND BERGER, “‘Encyclopedia of Physics,’’ Berlin, Springer-Verlag, Vol. 
38, Part 2, 1959. 
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Rapio WAVES IN THE IONOSPHERE, by K. G. 
Budden. 509 pages, 6 X 9in. New York, 
Cambridge University Press, 1961. Price, 
$18.50. 


The enormous success and expansion during 
the last twenty years of such fields as nuclear, 
quantum and solid state physics have kept 
many other important branches of science out 
of the limelight. The physics of the iono- 
sphere is one of the latter and, besides the 
electrical engineers, many physicists have 
maintained a great deal of interest and ac- 
tivity in this area. The public and inter- 
national interest in satellites and space travel 
has contributed to stimulate a new departure 
in this and other “space sciences.”’ 

At this juncture, Ratcliffe’s monograph on 
“The Magneto-Ionic Theory’? (Cambridge, 
1959) and the present book provide an up-to- 
date and comprehensive survey of this field 
of physics. While Ratcliffe’s works empha- 
size the physical aspects, the present treatise 
gives the full mathematical treatment of the 
problems involved, including those formulated 
in Ratcliffe’s. 

In the first half of the book the subject is 
introduced starting from the basic Maxwell 
and constitutive equations, and the propaga- 
tion of electromagnetic waves in homogeneous 
(isotropic and anisotropic) media. This first 
half deals with the phenomena where the 
“ray-approximation”’ is The second 
half gives the “full wave theory.”’ The sub- 
ject matter is developed carefully, all detini- 


valid. 


tions clearly stated and analyzed, standard 
definitions (for example, reflection and trans- 
The book is a 
mine of very diverse mathematical techniques 
(some, like the W.K.B. method, widely used 
in other fields) well exposed and adapted to 


mission coefficients) revised. 


the situation at hand. 

The research worker and the student alike 
will appreciate the methodical presentation, 
the careful bibliography, and the diversity 
of problems solved or proposed. 

HENRI AMAR 
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HANDBUCH DER Puysik (ENCYCLOPEDIA OF 
Puysics), VoL. X, STRUCTURE OF LIQUIDs, 
edited by S. Fliigge. 320 pages, diagrams, 

6} X 9} in. Berlin-Wilmersdorf, Springer- 

Verlag, 1960. Price, DM 76.80. 


While the physics and chemistry of solids 


have been experiencing an extraordinary 


growth and popularity, the investigations 
of the liquid state have been relatively limited 
in number and sometimes in scope. This is 
due partly to the fewer technological applica- 
tions of liquids and partly to the inherent 
difficulties of 


number — of 


and formidable the subject. 


Nevertheless a devoted and 
talented theoreticians and experimentalists 
have been pursuing an active research in this 


field. The 


stantial if not spectacular. 


success achieved has been sub- 

It is the purpose of this handbook volume 
to give a methodical survey of the present 
theories of the Structure of Liquids. This 
survey is divided into three parts. In Part 1, 
H. S. Green (well known for his early work 
with Max Born and for his book on the molec- 
ular theory of liquids) first analyzes the 
general nature of the liquid state, and the 
quantitative approach of the distribution 
function. 
on the structure of uniform, non-uniform and 
“quantum” Part II, 
Syn Ono (Radiation Laboratory, University 


of Tokyo) surveys the molecular theory of 


He then exposes the modern views 


liquids. In Professor 


surface tension in liquids. He describes 
clearly the three present approaches to this 
subject, namely, (a) through statistical me- 
through thermodynamics and 
“quasi-thermodynamics,” (c) through lattice 


The third and shortest part, by 


chanics, 


theories. 
F. P. Buff (University of 
voted to the Theory of Capillarity. 


Rochester), is de- 
Gibbs’ 
thermodynamic formulation is stressed and 
several applications are given. 

The nature of the subject makes it neces- 
sarily highly mathematical. The exposition 
is, on the whole, methodical and thorough, and 
the bibliography is excellent. The German- 


English dictionary appendix is hardly neces- 
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sary since the three articles are written in 
English. 
advantageously replaced by a comprehensive 


This appendix could have been 
and descriptive synthesis facilitating for the 
reader the long journey through the long 
mathematical treatments. 

Henri AMAR 
Temple University 


THE FERMI SURFACE, edited by W. A. Har- 
rison and M. B. Webb. 356 pages, dia- 
grams, 83 X 11 in. John Wiley & Sons, 
Inc., New York, 1961. Price, $10.00. 


This is the proceedings of an International 
Conference on the Fermi surface which was 
held at Cooperstown, N. Y., during August 
22-24, 1960. 


tended and the papers and discussions present 


The conference Was well at- 


a fair picture of the most recent thinking with 
regard to the electronic structure of metals 
Both 


are well represented in the proceedings. It 


alloys. theory and experiment 
is satisfying to find large and growing areas 


of agreement between the two. This is in 
marked contrast to the picture ten to fifteen 
years ago when it seemed that attempts to 
improve on the simple free electron theory of 
ended in agreement with 


metals poorer 


experiment. Today, theoretical methods of 
considerable power have given us some in- 
sight into why the free electron model is as 
good as it is. Another striking change is the 
abandonment of analytical forms for the shape 
of the Fermi surface and their replacement 
The 


the de Haas- 


by geometrical and topological models. 
major experimental techniques 
van Alphen effects, galvanomagnetic effects, 
cyclotron resonance, anomalous skin effect, 
magnetoacoustic effect, and transport effects 
—are described in separate sessions of from 
three to six papers each. Finally, the results 
of the conference are summarized by M. H. 
Cohen (theory) and A. B. Pippard (experi- 
ment). This is an excellent reference work 
whose usefulness is enhanced by the thought- 
fulness of the editors in providing a nine page 
subject index. 
F. J. DONAHOE 
The Franklin Institute Laboratories 
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FRAME ANALYsiIsS: A UNIFIED INTRODUCTION 
TO THE MATRIX ANALYSIS OF STRUCTURES, 
by A. S. Hall and R. W. Woodhead. 247 
New York, 

1961. Price, 


pages, diagrams, 6 X 9 in. 
John Wiley & Sons, 
$8.50. 


Inc., 


The electronic computer makes possible 
the analysis of complex framed structures by 
matrix methods. The authors of this concise 
yet comprehensive book successfully demon- 
strate that the representation in matrix form 
of a number of familiar methods of structural 
analysis leads to a more unified approach and 
helps to promote clearer understanding of 
the individual techniques. However, only 
an elementary knowledge of matrix methods 
is required, for which sufficient background 
is supplied in an Appendix. 

The presentation is divided into three parts. 
The first two introductory chapters summar- 
ize the basic definitions and concepts, such 
as equilibrium, compatibility, etc., and dis- 
tinguish as well as relate the two basic and 
reciprocal techniques for the analysis of in- 
determinate structures: Action methods (such 
as virtual work) in which restraints are re- 
placed by redundants of the proper magnitude 
to restore continuity; and displacement 
methods (such as slope deflection) which 
solve for the magnitude of the displacement 
required to restore equilibrium. Each of the 
two main parts of the book is based on one of 
these approaches. Their duality is demon- 
strated by solution of identical problems by 
both methods. 

The first two chapters of Part I discuss the 
calculation of influence and flexibility coeffi- 
cients and illustrate their use by a number of 
fully worked representative examples. Chap- 
ter 5 demonstrates the advantages of trans- 
forming the flexibility coefficients so as to 
produce more reliable numerical solution by 
change of scale factors, and of redundants; 
and Chapter 6 numerizes the mathematical 
technique for axis transformations. The final 
two chapters treat the flexibility analysis of 
a single member, and apply the techniques 
The 


use of the elastic center, and the column 


of transformation of restraints and axes. 


analogy, are shown to be particular cases of 
the general methods given. 

Part II, on Displacement Methods, is 
introduced by a discussion of the slope deflec- 
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tion method. The two succeeding chapters 
derive matrices for the stiffness of a single 
member and of a frame, and include a brief 
treatment of the effects of incomplete fixity. 
Chapter 12 reviews two relaxation methods 
related to slope deflection as applied to plane 
Moment “group 
which has been used where the 


frames: distribution, and 
relaxation,” 
presence of sidesway complicates moment 
distribution. <A final chapter deals with the 
method of coefficients of restraint, a form of 
moment distribution employing modified stiff- 
ness coefficients and carry-over factors which 
obviate the need for iteration. 

This book, written for the senior student 
and practicing engineer, is concerned with 
the application of mathematical method to 
structural analysis, and is not intended as a 
treatise on matrix calculus. The authors are 
to be highly commended for producing a very 
useful and instructive work within the limits 
which they have set themselves, and which 
they invite the reader to extend by further 
study of the references they supply 

K. L. CAPPEL 
The Franklin Institute Laboratories 


ADVANCES IN DOCUMENTATION AND LIBRARY 
ScriENCE, Vor. III, Part 1, INFORMATION 
RETRIEVAL AND MACHINE TRANSLATION, 
edited by Allen Kent. 686 pages, illustra- 
tions, 6 X 9 in. New York, Interscience 
Publishers, 1960. Price, $23.00. 


An estimated 70,000 scientific journals ap- 
pearing in many foreign languages will be 
published this year. Significant research re- 
sults reported in these publications will make 
vital contributions in every area of scientific 
and technical endeavors. One great problem 
that presently exists is that a large percentage 
of this information may never reach those 
The 
magnitude of the world’s printed scientific 


individuals who could make use of it. 


output coupled with the language barrier has 
created a communication problem that ap- 
pears will only be solved by new methods and 
techniques. These methods and techniques 
are directed towards the use of computers 
for the retrieval of information and the trans- 


lation of languages. 


REVIEWS 


Western Reserve University’s Center for 
Documentation and Communication Research 
has been one of the pioneering groups that has 
carried out research in this area. A series 
of volumes concerned with the many facets 
of this subject have been published over the 
past years. The present volume under review 
is based on papers presented at the Inter- 
national Conference for Standards on a Com- 
mon Language for Machine Searching and 
Translation, sponsored by Western Reserve 
University and Rand Development Corpora- 
tion held in Cleveland, Ohio, September 6-12, 
1959. The 21 
specialists drawn from all over the world, 


papers were presented by 


who were concerned with some methods of 
“working toward a common machine language 
or a series of compatible machine languages 
for processing scientific and technical litera- 
that it 
correlated, 


ture so may be searched, selected, 
and translated by automatic 
equipment.” The deal with 


1. Classification and coding of specitic 


papers many 
areas: 
types of literature for certain user groups, 2. 
Description of the specific computers and 
computer-like devices, 3. The problems of 
analysis of scientific information for the pur- 
pose of storage and retrieval and 4. Descrip- 
tion of specific systems that are now 
operational. 

This volume represents a significant review 
of the state of the art as it is now practiced in 
organizations located throughout the world. 
The principal criticism is that the papers re- 
flect a certain chaos that still exists in the 
information retrieval field. As Dr. Rangan- 
athan states in the appended discussion in 
this varies from 


volume, the terminology 


paper to paper with inconsistencies and lack 
The volume will not 


of general standards. 
be read from cover to cover by most readers 
but will stand as an important reference work. 
A more prudent editorial policy might have 
pared the size of the volume down to manage- 
able proportions and also deleted some of 
the over-lapping and duplication. Adequate 
author abstracts would have been of consider- 
able help in digesting some of the material. 
Joun P. McGowan 

The Franklin Institute Library 
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BOOK NOTES 


PHysicAL CHEMISTRY OF PROCESS METAL- 
LURGY (in two parts), edited by George 
R. St. Pierre. 1374 pages, illustrations, 
6X9 in. New York, Interscience Pub- 
lishers Inc., 1961. Price: $22.50 (Part 1); 
$25.00 (Part 2). 


These two volumes (Vols. 7 and 8 of the 
Metallurgical Society Conference series) are 
the proceedings of an international symposium 
sponsored by the AIME in 1959. Part 1 
deals with the physical chemistry of metal- 
lurgical and oxide phases, thermodynamics of 
metals, transport and mixing, the nature and 
structure of liquid metals, and solubility and 
phase equilibria in metals systems. 

art 2 covers process reaction rates and 
mechanisms, solidification of metals, proper- 
ties of halide and sulfide melts, industrial 
applications of principles, desulfurization in 
steel plants, and process control and statistical 
methods. 

Discussions are included with a majority of 
the individual papers. 


POLYMERIC MATERIALS, by Charles C. Wind- 
ing and Gordon D. Hiatt. 406 pages, dia- 
grams, illustrations, 6 X 9 in. New York, 
McGraw-Hill Book Co., Inc., 1961. 
$12.00. 


Price, 


The purpose of this book is to introduce 
the reader to polymeric materials so that he 
may be able to understand them and cooper- 
ate with others actively working in the fields 
of plastics, rubber, protective coatings and 
textile fibers. The book can also be used as 
a text for an introductory college course, but 
a good knowledge of chemistry is a_pre- 
requisite. It is not intended to aid the 
specialist in any of the polymeric fields 
covered. 


ENGINEERED CastTINGs, by Glenn J. Cook. 
240 pages, diagrams, illustrations, plates, 
6X9in. New York, McGraw-Hill Book 
Co., Inc., 1961. Price, $8.50. 


A practical approach to all phases of engi- 
neered castings is used in this new book aimed 


at design engineers, draftsmen, materials 
engineers, purchasing agents, and cost-reduc- 
tion specialists, as well as at foundrymen, pat- 
ternmakers, foundry engineers and managers. 
It covers how castings are designed, made 
and used, and gives practical suggestions for 
use in selecting cast material specifications 
and in buying castings. The author, Mana- 
ger of Marketing in GE’s Foundry Depart- 
ment, expresses the hope that the book will 
stimulate young men to find careers in the 
foundry industry. 


HEAT TRANSFER, by Benjamin Gebhart. 454 
pages, diagrams, 6 X 9 in. New York, 
McGraw-Hill Book Co., Inc., 1961. 
$10.75. 


Price, 


This book describes the more important 
physical processes, theories, and methods of 
analysis for the heat-transfer field. Although 
it is designed as a text, it will be of use to the 
practicing engineer as well as the student. 
Basic theories and formulations of conduction, 
radiation and convection are presented, and 
the mathematical statement of problems is 
emphasized. The author’s intent is to dem- 
onstrate the diverse and special considerations 
whereby general formulations are adapted to 
Laplace transforms 
are not used in the solution of partial dif- 


particular applications. 


For use as a text in a 
first course, the author suggests omission of 
much of Chapter 6 on fluid flow and all of 
Chapter 13 on experimental analogues. In 


ferential equations. 


some cases, the same symbol denotes several 
different quantities, which may prove a little 
confusing. 


Tue PuysicaL THEORY OF TRANSISTORS, by 
Leopoldo B. Valdes. 358 pages, diagrams, 
6} X 9} in. New York, McGraw-Hill Book 
Co., Inc., 1961. Price, $10.50. 


This book is intended for a graduate course 
in transistor physics, and as a reference for 
those already working in this expanding field. 
In order to unify the various concepts making 
up the study of transistors, some introductory 
material is included. 
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text analyzes the flow of current 


The 
through semiconductor materials and derives 
relationships between the electrical charac- 
teristics of transistors and their physical 
The book is divided into three 
parts, the first giving introductory or review 
The second part is devoted to the 


structure. 


material. 
more physical aspects of transistor theory 
and deals in generalities, and the third part 
analyzes specific device structures and the 


electrical characteristics of diodes and 
transistors. 
SPECIAL CERAMICS, edited by P. Popper. 363 


pages, diagrams, illustrations, 5} X 8} in. 

New York, Academic Press, Inc., 1961. 

Price, $10.50. 

This book represents the proceedings of a 
Symposium held by the British Ceramic Re- 
search Association in the summer of 1959. 
The symposium was held to give an oppor- 
tunity for information to be exchanged con- 
cerning the research on non-oxide ceramics 
conducted by the BCRA in cooperation with 
various government departments and indus- 
trial firms, in response to the growing demand 
for new high temperature materials. 
dealt with at the Symposium included: Prop- 
erties and Structure, Measurements; Prepara- 
tion and Properties of Nitrides; Preparation 
and Properties of Other Non-Oxides; 
Furnaces, Techniques, Analyses, etc. 


Topics 


BIOLOGICAL 
1119 
pages, diagrams, illustrations, 7 X 10 in. 
New York, Reinhold Publishing Corp., 
1961. Price, $20.00. 
The biological sciences (from Abiogenesis 
to Zoogeography) are covered in over 800 
articles specially written for this new encylo- 


THE ENCYCLOPEDIA OF THE 
ScIENCEs, edited by Peter Gray. 


pedia. Its up-to-dateness is indicated by 
the fact that it includes articles on biological 
warfare (by LeRoy Fothergill) and on space 
biology (by O. E. Reynolds and John M. 
Talbot). It deals with developmental, eco- 
logical, functional, structural 
taxonometric aspects of the biological sciences, 


genetic, and 
and covers several topics in biophysics and 
The contributors include ex- 
References are 
given at the end of most of the articles. The 
expert and the layman alike will find answers 
to their questions in this newest of the special- 


biochemistry. 
perts from many countries. 


ized encyclopedias. 


PowpDER METALLURGY, edited by Werner 
Leszynski. 843 pages, diagrams, illustra- 
tions, 6 X 9 in. New York, Interscience 
Publishers, 1961. Price, $25.00. 

The papers in these proceedings were pre- 

Powder Metal- 

Four main 


sented at the International 
lurgy Conference in June, 1960. 
divisions cover: Mechanism of Sintering (13 
papers) ; Dispersion Strengthening (6 papers) ; 
Technology and Methods (12 papers); and 
Materials and Applications (12 papers). 
Discussions are included with 29 of the papers. 
Paul Schwarzkopf (one of the pioneers in 
this field) gave the keynote address, which is 
included as an introduction to the volume. 
It covers the last fifty years of the powder 
metallurgy art. The papers are about equally 
divided between theoretical and practical ap- 
proaches, and deal with powder metallurgy 
as a science, a technology, and an industry. 
The stature of the contributors and the wide 
range of material covered make this a fine 
and lasting contribution to the literature in 
this important field. 


QuantuM MECHANICs, by Eugen Merzbacher. 

544 pages, diagrams, 6 X 9in. New York, 
John Wiley & Sons, Inc., 1961. Price, 
$12.00. 


With emphasis on contemporary develop- 
ments, this book represents a new and valu- 
able contribution to the subject of non-rela- 
tivistic quantum mechanics. The develop- 
ment of the discussion of quantum mechanics 
is presented in three stages by the author. 
In the first, the ordinary wave mechanics is 
developed, enabling the student to handle 
actual problems of atomic and nuclear physics 
as early as possible. The matrix form of 
quantum mechanics is explicitly introduced 
as a useful tool in formulating one dimensional 
problems and is used again in the treatment of 
spin, in the second stage. In the third stage, 
the two versions of quantum mechanics are 
unified in the general and more abstract 
Dirac formulation which the author employs 
throughout the second half of the book 


FUNDAMENTALS OF MOopeERN Puysics, by 
Robert M. Eisberg. 729 pages, diagrams, 
6 X9in. New York, John Wiley & Sons, 
Inc., 1961. Price, $10.50. 


This text provides a rapid but thorough 
treatment of the important developments 
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. 
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leading up to Schroedinger quantum me- 
chanics; a careful introduction to the basic 
theory of that subject, using an inductive 
rather than a postulatory approach; and il- 
lustrations of the theory through a complete 
discussion of a large number of typical calcu- 
lations. Finally, he uses the theory to give 
a detailed explanation of the properties of 
atoms and nuclei. The author launches each 
topic near its beginning by means of arguments 
that presuppose only a background in ele- 
mentary physics, and in mathematics through 
intermediate calculus. 


ENGINEERING DATA FOR PropuCT DESIGN, 
edited by Douglas C. Greenwood. 425 
pages, diagrams, 8} 11} in. New York, 
McGraw-Hill Book Co., Inc., 1961. Price, 
$10.00. 


This manual has been prepared to help the 
engineer solve design problems. Through 
charts, formulas, tables, and special articles 
(taken from Product Engineering) the author 
places valuable and time-saving data at the 
fingertips of designers, about metals and al- 
loys, nonmetallic materials, finishes, and de- 
sign analysis. Clearly illustrated and care- 
fully indexed, the manual should be of real 
help in answering questions about beams, 
torque, gears, shafts, control, and the appli- 
cation of electricity, electronics, hydraulics, 
and pneumatics to product engineering. 


ELECTRICAL ENGINEERING, by S. B. Ham- 
mond. 547 pages, diagrams, 7 X 9} in. 
New York, McGraw-Hill Book Co., Inc., 
1961. Price, $8.75. 


This text is intended to be used for engi- 
neering students not majoring in electrical 
It covers the basic principles of 
Sample chap- 


engineering. 
electricity, its use and control. 
ter headings are: Characteristics of Voltages 
and Currents; Electric Networks; Nonlinear 
Resistor Circuits; Electron Physics ; Feedback 
Circuits; Linear Approximations to Machine 
Analysis; Digital Computers; and Magnetic 
Fields and Circuits. The three main divisions 
deal with electric circuit analysis, electronics, 


and conversion of energy. 


Book Notes 
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DESIGN AND CONSTRUCTION OF PORTS AND 
MARINE StructTuRES, by Alonzo DeF. 
Quinn. 521 pages, illustrations, 6 X 9 in. 
New York, McGraw-Hill Book Co., Inc., 
1961. Price, $16.00. 


In this work, the author brings together 
the essential concepts of designing and build- 
ing modern shipping-terminal, harbor, and 
offshore marine structures—along with the 
working principles and practices that can be 
adapted to specific project requirements. He 
provides you with authoritative construction 
know-how and recom mendations for achieving 
top efficiency and economy in meeting the de- 
mands of today’s advanced engineering tech- 
This is the only U. S. published book 


niques. 

dealing solely with ports and marine 
structures. 

FUNDAMENTALS OF CiRCUIT THEORY, by 


Norman Balabanian. 548 pages, diagrams, 
6 X 9 in. Boston, Allyn and Bacon, Inc., 
1961. Price, $9.50. 


The author, Norman Balabanian of the 
Electrical Engineering Department of Syra- 
cuse University, states in his preface that this 
book is his attempt at giving an account of 
the concepts and techniques of circuit theory. 
The book’s 13 chapters include Resistance 
Networks, Complete Response of Simple Net- 
works, Sinusoidal Steady-State Response, 
The Complex-Frequency Variable, Two-Port 
Networks, Electronic and Their 
Models, Operational Diagrams and Signal 
Flow Graphs, and Complete Response to 
Arbitrary Wave Forms. 

Intended primarily as a the 
author has attempted to leave the student 
with the impression that circuit theory can 
still challenge his inventiveness, that it 
still has plenty of vitality and that it is not, 
by any means, a finished subject. 


Devices 


textbook, 


AN INTRODUCTION TO THE THEORY AND 
PRACTICE OF ‘TRANSISTORS, by F. F. 
Roberts and J. R. Tillman. 335 pages, 
diagrams, plates, 53 X 8} in. New York, 
John Wiley & Sons, Inc., 1961. 
$8.00. 


Price, 


The authors have prepared this book to 
present a wide treatment of the subject in 
one volume, for newcomers to the transistor 
and related semiconductor 


fields. Part 1, 
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154 Book 
covering basic physical theory, will be inter- 
esting to electronics engineers and experi- 
mental physicists; Part 2 covers technology, 
properties and applications. 
sistors are omitted, as are the merits of dif- 


Unipolar tran- 


ferent types of sealing and the factors af- 
fecting reliability, although a brief section is 
devoted to the point-contact transistor. 
MaAGNetic CONTROL OF INDUSTRIAL Morons, 

by Gerhart H. Heumann. Part 1: A-C 
ContTROL Devices AND ASSEMBLIES, 273 
pages, illustrations, 6 X 9 in. New York, 
John Wiley & Sons, Inc., 1961. 
$9.00. Part 2: A-C Motor CONTROLLERS, 
325 pages, illustrations, 6 X 9 in. New 
York, John Wiley & Sons, 1961. 
Price, $9.00. 


Price, 


Inc., 


This three-volume work divides the subject 


of magnetic control of industrial motors 
logically into three parts, two of which are 
reviewed here. The third D-C Motor 
Controllers) was not available for comment 


at the time this issue of the JOURNAL went to 


(on 


press. Primarily an application book, the 
work covers magnetic devices and accessories 
used for control of industrial motors. It 
deals with principal circuits, selection of con- 
troller sizes and components, economic factors 
affecting controller selection, motor protec- 
tion, and safety codes and standards. 

Part 1 includes chapters on: Symbols and 
Diagrams; Alternating-Current Contractors 
and Relays; Accessories; Control-Circuit De- 
vices; Resistors Rheostats; Solenoids, 
Thrustor Mechanisms, and Brakes; Mechani- 
cal Design of Controllers; and Maintenance. 

Part 2 covers: Drive Characteristics ; Con- 
trol of Squirrel-Cage Motors; Control of 
Wound-Rotor Motors; Logic Functions and 
Static Switching; Motors Protection; 
Motor Branch-Circuit Protection. 


and 


and 


INTRODUCTION TO ComPLEX ANALYsiIs, by 
Zeev Nehari. 254 pages, diagrams, 6 XK 9 
in. Boston, Allyn & Bacon, Inc., 1961. 
Price, $7.50. 


This text provides an introduction to the 
basic techniques and methods of complex 
analysis. Since these techniques may be 
applied in many areas, the book should be 
students in mathematics, 
Suitable for a 


useful for 
engineering and the sciences. 


many 


Nores J. 


senior or first-year graduate course, the ma- 
terial presented covers complex variables, 
analytic functions, complex integration, the 
application of Cauchy’s theorem, conformal 
mapping, physical applications such as bound- 
ary-value problems and two-dimensional fluid 
flow. Sections on analytic continuation, the 
symmetry principle and the harmonic conju- 
gate may be omitted if it is desired to place 


emphasis on applications. 


SEMICONDUCTOR DEVICES AND APPLICATIONS, 
A. Greiner. 487 pages, diagrams, 
6} X in. New York, McGraw-Hill Book 
Co., Inc., 1961. Price, $12.50. 


This new textbook, suitable for the senior 
or starting graduate student, emphasizes the 
proposition that an engineer must have a 
deeper understanding of the operation of 
electron devices than merely the ability to use 
The book con- 


by R. 


specific devices in circuits. 
tains a thorough discussion of the properties 
of semiconductors and junctions as a prelimi- 
nary to the discussion of transistor operation. 
Characteristics and parameters are derived 
from basic physical concepts, and extensive 
application of transistors to a variety of small 
and large signal circuits is included. The 
author has inserted many completely worked 
examples. 


INTRODUCTORY SYSTEM ANALYsIS: SIGNALS 
AND SYSTEMS IN ELECTRICAL ENGINEERING, 
by William A. Lynch and John G. Truxal. 
445 pages, diagrams, illustrations, 6 XK 9 
in. New York, McGraw-Hill Book Co., 
Inc., 1961. Price, $7.50. 

The primary objective of this new text is 
to present the and 
techniques of modern electrical engineering 
by introducing the viewpoints, approaches, 
and methods of the electrical engineer, and to 
familiarize the reader with electrical devices 


fundamental concepts 


and equipment. It also provides the requi- 
site foundation for future study, and presents 
a definitive picture of the role of the electrical 
engineer in the modern interdisciplinary sys- 
tems engineering efforts, both military and 
industrial. 

In addition to the basic techniques of sys- 
tem analysis being presented simply, but 
thoroughly, there are also several unusual 
topics treated, among them signal-flow graphs 


and analog computers. 
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PUBLICATIONS RECEIVED 


RADIOACTIVE SUBSTANCES, by M. Curie. A translation from the French of the classical thesis 
presented to the Faculty of Sciences in Paris. 94 pages, 44 X 7} in. New York, Philoso- 
phical Library, 1961. Price, $2.75. 


THt WorLb oF GEOLOGY, edited by L. Don Leet and Florence Leet. 262 pages, illustrations, 
53 X 8$in. New York, McGraw-Hill Book Co., Inc., 1961. Price, $4.25. 


A MIpWESTERN ALMANAC, by Austin L. Rand and Rheua M. Rand. 176 pages, diagrams, 
6 X9in. New York, The Ronald Press Company, 1961. Price, $4.50. 


TEXT-BOOK ON SPHERICAL AsTRONOMY, by W. M. Smart. Fourth edition, reprinted, 430 
pages, diagrams, 5$ K 8} in. New York, Cambridge University Press, 1961. Price, 
$3.95 (paper). 


ELEMENTS OF MATHEMATICS, by J. Houston Banks. Second edition, 465 pages, diagrams, 
6 X 9in. Boston, Allyn and Bacon, Inc., 1961. Price, $6.95. 


HuMAN HerREpITY, by Jean Rostand. Translated from the French by Wade Baskin. 135 
pages, 44 X 7{ in. New York, Philosophical Library, 1961. Price, $4.75. 


Second revised edition, 1333 pages, plates, 
Price, 84s. net. 


PuysicAL CHEMIstTrRY, by E. A. Moelwyn-Hughes. 
diagrams, 6 X 9} in. London, Pergamon Press, 1961. 


MicwHAEL FARADAY. A List oF His LECTURES AND PUBLISHED WRITINGS, compiled by Alan 
E. Jeffreys, with a foreword by Sir Lawrence Bragg. 86 pages, plates, 7 X 92 in. New 
York, Academic Press Inc., 1961. Price, $7.00. 


Third edition, 437 pages, plates, diagrams, 7} X 10 
Price, $12.50. 


ELECTRICAL EsTIMATING, by Ray Ashley. 
in. New York, McGraw-Hill Book Co., Inc., 1961. 


ELEMENTARY FLuIp MECHANICS, by John K. Vennard. Fourth edition, 570 pages, diagrams, 
6 X9in. New York, John Wiley & Sons, Inc., 1961. Price, $7.95. 


INTRODUCTION TO NUCLEAR ENGINEERING, by Raymond L. Murray. Second edition, 384 
pages, illustrations, diagrams, 6 X 9 in. Englewood Cliffs, Prentice-Hall, Inc., 1961. 
Price, $12.00 (trade) ; $9.00 (classroom adoption). 


188 pages, diagrams, illustrations, 53 X 8} in. 
Price, $3.20 (paper) ; $5.00 (cloth). 


INTRODUCTION TO H1-F1, by Clement Brown. 
New York, Gernsback Library, Inc., 1961. 


PROBLEMS IN QuANTUM MEcHANICS, by I. I. Gol’dman, V. D. Krivchenkov, V. I. Kogan, 
and V. M. Galitskii, translated, edited and arranged by D. ter Haar. 394 pages, diagrams, 
6 X 9in. New York, Academic Press Inc., 1961. Price, $12.00. 


TELEVISION AND RApIto REPAIRING, by John Markus. Second edition, revised and up-dated, 
568 pages, illustrations, 6 X 9 in. New York, McGraw-Hill Book Co., Inc., 1961. 


Price, $8.95. 


CyBERNETICS, by Norbert Wiener. Second edition, revised, 212 pages,6 X 9 in. New York, 
John Wiley & Sons, Inc., 1961. Price, $6.50. 
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NEW BOOKS IN THE FRANKLIN INSTITUTE LIBRARY 


PHYSICS 


FirtLD THEORY FOR ENGINEERS, by P. Moon. Van Nostrand, 1961. 

CompuTING METHODS AND THE PHASE PROBLEM IN X-RAY CRysTAL ANALYsIS, edited by 
Ray Pepinsky. Pergamon, 1961. 

THE ENCYCLOPEDIA OF Microscopy, by G. L. Clark. Reinhold, 1961. 

AN INTRODUCTION TO THE THEORY OF VIBRATING SysTEMS, by W. G. Bickley. Clarendon, 
1961. 

MATERIALS FOR NUCLEAR ENGINEERS, edited by A. B. McIntosh. Interscience, 1960, 

Puysics OF SEMICONDUCTORS, by A. F. Ioffe. Academic, 1960. 

PHYSICS FOR THE INQUIRING MIND, by E. M. Rogers. Princeton University Press, 1960. 

GENERAL RELATIVITY AND GRAVITATIONAL WAVES, by J]. Weber. Interscience, 1961. 

TECHNIQUES OF HIGH ENERGY Puysics, by D. M. Ritson. Interscience, 1961. 

X-Ray Microscopy, by V. E. Cosslett. Cambridge University Press, 1960. 

THE ROTATION OF THE EARTH, by W. H. Munk. Cambridge University Press, 1960. 

DISPERSAL RELATIONS AND CAUSAL DescRIPTION, by J. Hilgevoord. North-Holland, 1960. 

MAGNETO-FLUID Dynamics, edited by F. N. Frenkiel. National Research Council, 1960. 

QuantuM MEcuHANICs, by J. L. Powell. Addison-Wesley, 1961. 


CHEMISTRY AND CHEMICAL ENGINEERING 


CERAMICS, by P. W. Lee. Reinhold, 1961. 

PORCELAIN ENAMELS, by A. I. Andrews. 2nd ed. Garrard Press, 1961. 

AN INTRODUCTION TO TRANSITION-METAL CHEMISTRY: LIGAND-FIELD THEORY, by L. E. Orgel. 
Wiley, 1960. 

INTRODUCTION TO PETROLEUM CHEMICALS, by H. Steiner. Pergamon, 1961. 

PuysIcAL CHEMISTRY, by W. F. Sheehan. Allyn and Bacon, 1961. 


CIVIL ENGINEERING 


AMERICAN RAILROADS, by J. F. Stover. University of Chicago Press, 1961. 
ADVANCED DESIGN IN STRUCTURAL STEEL, by J. E. Lothers. Prentice-Hall, 1960. 

THE DESIGN OF CYLINDRICAL SHELL Roors, by J. E. Gibson. 2nd ed. Van Nostrand, 1961. 
PAPERS ON BUILDING ConstRUCTIONS. American Society for Testing Materials, 1960. 
SYMPOSIUM ON TECHNICAL DEVELOPMENTS IN THE HANDLING AND UTILIZATION OF WATER AND 
INDUSTRIAL WASTE WATER. American Society for Testing Materials, 1960. 


ELECTRICAL ENGINEERING 


LINEAR GRAPHS AND ELECTRICAL Networks, by S. Seshu. Addison-Wesley, 1961. 

TRANSISTOR Circuit ANALysis, by M. V. Joyce. Addison-Wesley, 1961. 

ESSENTIALS OF DIELECTROMAGNETIC ENGINEERING, by H. M. Schlicke. Wiley, 1961. 

ELECTRONICS, by P. M. Chirlian. McGraw-Hill, 1961. 

SEQUENTIAL DECODING, by J. M. Wozencraft. M.I.T. Press, 1961. 

ELSEVIER’S DICTIONARY OF AUTOMATION, COMPUTERS, CONTROL, AND MEASURING IN SIX 
LANGUAGES. Elsevier, 1961. 

ADAPTIVE CONTROL SysTEMS, edited by E. Mishkin. McGraw-Hill, 1961. 

Rapio WAVES IN THE IONOSPHERE, by K. G. Budden. Cambridge University Press, 1961. 
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LECTURES ON COMMUNICATION SysTEM THEORY, by E. J. Baghdady. 


McGraw-Hill, 1961. 
TRANSMISSION OF INFORMATION, by R. M. Fano. M.I.T. Press, 1961. 


METALLURGY 
X-Ray METALLOGRAPHY, by A. Taylor. Wiley, 1961. 
THEORY OF MetaL CurttinG, by P. H. Black. McGraw-Hill, 1961. 
RARE METALS HANDBOOK, by C. A. Hampel. 2nd ed. Reinhold, 1961. 
TECHNOLOGY OF NON-METALLIC COATING, by A. Ya. Drinberg. 
COLLOQUE DE ME&TALLURGIE: CORROSION, SECHE ET AQUEUSE. 

de Saclay, 1960. 

METALL UNTER DEM MIKROsKOP, von E. Kauczor. Springer, 1960. 


Pergamon, 1960. 
Centre d’études nucléaires 
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CURRENT TOPICS 


New Superconductor Material.—.\ 
new superconductor material which 
offers no resistance to electrical cur- 
rent has been developed by Atomics 
International, a division of- North 
American Aviation, Inc., under a re- 
search contract with the Atomic En- 
ergy Commission. 

A wire which can be easily coiled 
has been drawn from the new super- 
conductor. At liquid helium tempera- 
tures (—452 degrees F.) the wire has 
conducted 100,000 amperes of electri- 
cal current per square centimeter in 
a moderately high magnetic field 
(30,000 gauss). The performance of 
this wire may be roughly equated with 
operation of a household electrical ap- 
pliance at about 100,000 watts from 
a normal-size (110 volt) wall outlet. 

The new material is a cold-worked 
alloy made up of approximately three 
parts of niobium to one part of zir- 
conium. The alloy is malleable and 
strong and can be made into wires, 
bars, strips and other shapes without 
losing its superconducting properties. 
It is expected to retain these charac- 
teristics in substantially higher mag- 
netic fields. Its development makes 
possible the fabrication of supercon- 
ductor magnets to replace conven- 
tional bulky and relatively inefficient 
iron-core electromagnets. The super- 
conductor magnet has no iron core and 
may produce magnetic fields many 
times higher than conventional elec- 
tromagnets in common application 
today. 

Potentially, the use of supercon- 
ducting magnets of 100,000 gauss 
strength for high energy particle ac- 
celerators, or in connection with con- 
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trolled thermonuclear research de- 
vices, could substantially reduce their 
cost of operation. Similarly, super- 
conducting magnets may be particu- 
larly important in space technology 
where light-weight magnets are 
needed. Such magnets could store 
more electrical energy in less space 
than conventional electrical capacitors. 


New Superpressure Balloon De- 
veloped.— The development and suc- 
cessful testing of a new type of re- 
search balloon has been announced by 
the Air Force Cambridge Research 
Laboratories. The most significant 
feature of the development is that 
the balloon maintains a constant alti- 
tude without expenditure of ballast. 

Key to the design of the balloon is 
that its internal pressure exceeds that 
of the surrounding air. Research on 
such balloons, called superpressure 
balloons, has been under way for 
several years by scientists at AFCRL’s 
Geophysics Research Directorate. 

Tests of the newly designed balloon 
in April were the most successful yet 
conducted with superpressure bal- 
loons. The flight of longest duration 
took place on April 17. The balloon 
was spherically shaped with a diam- 
eter of 34 ft. Construction was of a 
Mylar material one and a half mils 
thick (one and a half thousands of an 
inch). 

The balloon was launched from 
AFCRL’s Vernalis, Calif. test site. 
It carried a payload of 40 Ib. and was 
designed to float at an altitude of 
70,000 ft. Nine days later the test 
was completed and the flight auto- 
matically terminated over the ocean 
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400 miles off the coast of the State of 
Washington. 

During the nine-day flight, the bal- 
loon maintained a constant volume 
and consequently a constant altitude. 
This was achieved by inflating the 
balloon with additional gas to create 
an overpressure condition at floating 
altitude. At night the pressure in 
the balloon decreased due to the loss 
of radiative heating from the sun, but 
the volume of the non-extensible bal- 
loon remained the same. Conven- 
tional zero-pressure balloons decrease 
in volume after sunset and require the 
release of ballast to maintain floating 
during the night. In the case of long 
duration flights the amount of ballast 
required can exceed the gross weight 
of the balloon and its payload. 

The flight duration was a record for 
any type of superpressure balloon. 
The actual floating altitude was 70,500 
ft. Data telemetered from the balloon 
throughout the flight indicate that the 
balloon could have continued its flight 
for a period considerably longer than 
the nine days. 

The results of the test series prove 
the feasibility of superpressure bal- 
loons of this size and payload capacity. 
A significant step has been achieved 
toward the development of higher fly- 
ing, larger payload capacity super- 
pressure balloons that will carry 100 
pounds to a constant level altitude of 
100,000 ft. for periods of 10 days or 
more. 


Camera Photographs Planets, Satel- 
lites in Daylight.— A special Air Force 
camera that looks like an over-size 
anti-aircraft gun has been photograph- 
ing missiles, satellites and planets in 
broad daylight. 

Consisting of 19 long-barrelled 5-in. 
refracting telescopes linked to as many 
TV-like image-orthicon tubes, the 
giant “‘Facet-Eye Camera” has been 
installed at the Air Force Missile De- 
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velopment Center, Holloman Air Force 
Base, New Mexico. The center is one 
of seven maintained by the Air Force 
Systems Command. 

Moving on a massive but delicately- 
balanced tracking mount, the system 
recently captured continuous bright- 
image views of Venus and Jupiter in 
broad daylight. 

Planetary images—not visible to 
the naked eye—were transmitted to 
indoor viewing screens where they 
were photographed with a precision- 
plate camera to produce clear, distor- 
tion-free pictures ordinarily obtained 
only during night hours. The total 
field covered is many times the size 
obtainable standard tracking 
telescopes. 

In terms of daylight star detection, 
this represents an advance of approxi- 
mately six stellar magnitudes. Scan- 
ning the heavens at night, the new 
camera system has photographed a 
star cluster, The Pleiades, with image 
resolution of faint stars up to the 12th 
magnitude. 

The number of telescopes and oscil- 
loscope screens will be increased to 
25 to complete the planned design of 
the Facet-Eye and to provide a view- 
ing field of 2} by 7 degrees. 

Basically conceived for the purpose 
of tracking missiles, satellites and 
other space-bound objects at extreme 
distances and in any kind of light, 
the system developed in 
AFMDC’s orbital Mechanics Divi- 
sion by Walter E. Woehl, theoretical 
and applied physicist. Woehl for- 
merly held high research positions in 
his native Germany and at the Royal 
Arsenal in England. 

Described as the only instrument of 
its kind in the world, the system was 
built to exacting specifications by the 
International Telephone and Tele- 
graph Laboratories, Fort Wayne, Ind. 
The 25 telescopes have been furnished 
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by the E. C. Larr Optics Company, 
Altadena, Calif. The 13,000-Ib., four- 
axis tracking mount built by 
Joseph Nunn and Associates, Pasa- 


dena, Calif. 
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New Information System. Comb- 
like arms flying on layers of air can 
monitor as many as 280 million char- 
acters of data in an information stor- 
age system recently announced by 
International Business Machines 
Corporation. 

The new 1301 disk storage unit can 
be linked to any of five of IBM's inter- 
mediate to large solid-state computers 
or shared by any two of them. It 
transforms them into advanced infor- 
mation centers capable of relaying 
constantly updated information, by 
means of IBM 
vices, to remote 
sandths of a second. 

Each 1301 unit holds from 50 to 56 
million characters of information de- 
pending on the computer employed. 
A total of five units can be used with 
any of the five computers. 

“The 1301, with its vast capacity 
and very high access speeds, opens 
the way to total systems processing,’ 
said Warren C. Hume, president of 
IBM’s Data Processing Division. 

“This means the vital records of 
most businesses can be maintained at 
a single, convenient location, provid- 
ing anyone within the company 
whether at home office, branch office 
or plant—almost immediate access to 
all pertinent data.” 

“New methods of data organiza- 
tion will enable our major computers 
to process records in any order without 
being slowed while files are searched. 
This significantly extends the range 
of IBM’s in-line random access equip- 
ment and makes possible many new 
processing concepts.”’ 

The 1301 is available for use with 
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the IBM 1410, 7070, 7074, 7080 and 
7090 computers. The file may be 
shared by any two of these systems. 

Each file contains one or two mod- 
ules (or stacks) of 20 rapidly revolving 
disks for data storage. ‘The 1301 has 
a read/write head for each disk sur- 
the end of 
comb-like access arms, are aligned in 
a vertical cylinder. Thus, at any 
given moment, they are positioned op- 
posite corresponding tracks or grooves 
on all forty disk surfaces in a module. 

This arrangement makes possible 
the picking up and depositing of data 
in corresponding data tracks with no 
access movement. Related informa- 
tion needed for frequent reference can 
parallel tracks 
forming concentric cylinders of data 
very fast reference. The only 
motion required within the file’ is 
horizontal from data cylinder to data 
cylinder. There is no necessity for 
vertical motion by the access arms, 
as in previous files, because there is 
a head for every disk surface. 

While reading in the same data 
cylinder, as many as 102,000 charac- 
ters in a module, 204,000 characters 
in one disk storage unit, or 1,020,000 
characters in a maximum system of 
five units may be reached with no 
access motion. These figures vary 
slightly with the computer used. The 
access mechanism requires a minimum 
of fifty milliseconds (thousandths of 
a second) to move horizontally from 
one cylinder to the next cylinder 
within small groupings of tracks on 
a disk surface. Within larger group- 
ings of fifty tracks, maximum access 
time is 120 milliseconds. The maxi- 
mum access motion time to any data 
in the 1301 is 180 milliseconds. Rec- 
ord organization techniques can cut 
this time substantially. 

The rate at which the 1301 reads 
characters in sequence into core stor- 
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age is between 75,000 and 90,000 
characters a second, depending on the 
type of system to which it is linked. 
This is faster than many tape systems. 
It is an important factor in sequential 
jobs such as recalling tables or pro- 
grams stored on the disks into core 
memory ; or transferring of the file to 
magnetic tape for a more permanent 
record of the data. 

A flexible record length feature adds 
to the efficiency of the 1301. Records 
(units of information) of different 
lengths may be stored on the cylinders. 
The capacity of the file, therefore, is 
much greater than that of fixed length 
files in which records are placed in 
identical storage spaces regardless of 
size. 

The read/write heads are part of a 
gliding shoe designed to fly along the 
disk surface. These shoes actually 
glide on a layer of air created by rota- 
tion of the disk at 1,800 rpm. There 
is an equilibrium of forces which keeps 
the head balanced less than a hair’s 
breadth above the disk surface. 

The 1301 makes feasible many 
large-size or large-volume jobs not 
previously possible within the com- 
puting range of the associated systems. 

Such giant tasks as keeping track 
of the inventory of a large mail order 
house or the updating of the employ- 
ment security accounts kept by the 
state government are within the scope 
of the 1301 files. 

Typically, the 1301 can be used as 
the basis for an information system 
for insurance. The detailed policy 
records, loan accounts and dividend 
accounts of a life insurance company 
can be stored in the file. The data 
would serve the computer for the 
processing of billing, anniversary cal- 
culations and renewals. In addition, 
it would also be immediately available 
anywhere in the company for policy- 
holder and agent service. 
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The loan department, for example, 
could ask through a remote unit a 
half-mile away from the central com- 
puter the immediate status on policy 
equities for loan approval on any one 
of thousands of policies. 

As well as extending the jobs a 
computer can do, the 1301 is capable 
of increasing the over-all efficiency of 
present systems. 

All of the programs, or job instruc- 
tions, for an IBM 7090 can be stored 
on a 1301 unit. Instead of the usual 
several minutes required to set up a 
7090 for a two or three minute proces- 
sing job, the 1301 can rush another 
program into the 7090 in seconds, 
virtually closing time-consuming and 
costly gaps between jobs. Sharing 
the same 1301 file, the scientific 7090 
and commercial 7080 can have refer- 
ence to a wide range of both scientific 
and commercial programs. 

Developed and manufactured at 
IBM facilities at San Jose, California, 
the 1301 will be shipped to customers 
starting in the third quarter of 1962. 


GaAs Detector Provides Greater 
Sensitivity.—The Lansdale Division 
of Philco Corporation has announced 
a gallium arsenide photovoltaic de- 
tector in a T0-18 case which marks 
the first exploitation of the small sig- 
nal photo detection capabilities of a 
solid state device. 

Reflecting the growing emphasis on 
miniaturization and reliability of elec- 
tronic components used in space sys- 
tems, Philco has placed the new single 
crystal photovoltaic detector in the 
small, rugged T0-18 case, used com- 
monly as a transistor package. 

Developed by engineers at the Lans- 
dale Division’s Special Products op- 
eration, the new device, GAU-401, 
provides high sensitivity for both 
visible and near-infrared radiation de- 
tection for the first time. 
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The detector is designed for the 0.4 
to 0.9 micron region, having a peak 
sensitivity at 0.85 microns. 

Dr. C. H. Sutcliffe, general manager 
of the special products operation, ex- 
plained that the new device provides 
greater sensitivity than do related 
photodevices such as photomultiplier 
tubes and silicon photodiodes. 

“This detector,” he said, “unlike 
conventional photo diodes, does not 
suffer serious sensitivity degradation 
when exposed to daylight and_re- 
places photo multipliers in providing 
1 to 2 orders of magnitude of increased 
responsivity in the near infrared 
region.” 

Dr. Sutcliffe emphasized the de- 
vice’s ability to sense during daylight 
as being most significant to celestial 
navigation of satellites and space 
vehicles. 

He explained that due to a high 
energy gap the GAU-401 can exhibit 
very high sensitivity at even higher 
temperatures than those at which 
silicon devices operate. Operation up 
to 120° C has shown no loss in de- 
tector sensitivity, he said. 

He also indicated that the device 
exhibits virtually no microphonism 
and can withstand massive radiative 
overloads without permanent damage. 

Dr. Sutcliffe stated that optimum 
device characteristics are observed at 
room temperature, obviating the need 
for dewar-cryostat arrangements asso- 
ciated with conventional IR detectors. 

Philco’s use of a modified TO0-18 
case (used for miniature transistors) 
marks an innovation in detector pack- 
aging, providing size reduction and 
ruggedness. 

It is estimated that the GAU-401 
detector is 1/20 the size of a related 
photomultiplier. 

Dr. Sutcliffe explained that the 
miniaturization of the GAU-401 de- 
vice promises to make possible the 
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combination of a solid state amplifier 
and photo detector potted in one single 
block no larger than a matchbox. 

The conventional T0-18 case has 
been modified with a glass window 
sealed at the top of the cylindrical 
package which measures 0.18 inch 
in diameter and 0.20 inch in length 
(excluding leads). 

In the standard configuration, 
square or circular cell areas ranging 
from 0.2 to 4.0 mm? can be provided. 

Dr. Sutcliffe stated that the GAU- 
401 becomes the first visible-near in- 
frared detector specifically designed 
for small signal use. 

He added that the wide variety of 
applications which the new detector 
is suited to include: celestial body 
sensing, missile tracking, fuzing and 
warning, space vehicle and _ satellite 
attitude control and guidance, stellar 
supervision of inertial platform, com- 
puter readout and numerous industrial 
process controls. 


Sapphires to Protect Phone Satel- 
lites.—Thousands of pieces of man- 
made sapphire will cover the surface 
of the communications satellites now 
being developed by Bell Telephone 
Laboratories, according to the com- 
pany’s president, Dr. James B. Fisk. 

The sapphires are expected to pro- 
tect solar cells from space radiation, 
enabling ‘“‘working’’ telephone satel- 
lites to endure the rigors of space for 
10 years or more. 

Dr. Fisk cited the sapphire protec- 
tion as an example of thorough-going 
measures that are being taken to 
achieve long-life satellite reliability. 
Bell Laboratories is using an approach 
it evolved in developing submarine 
cable systems. Some 1500 electron 
tubes have functioned under the 
oceans for the past two to six years 
without a failure. 

“Now the exacting methods that 
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produced extremely reliable submarine 
cable systems are being applied to 
satellites,” he said. ‘‘The result 
should be a shipshape, spaceworthy 
communications relay station.” 

He pointed out, however, the pres- 
sing need for a “‘real-life test’’ in space. 

Until now, most satellites have been 
designed for short-term experiments in 
space. Bell engineers have set as their 
goal the production of satellites that 
will work dependably for many years, 
serving as radio relay stations for 
communications overseas—telephone 
calls, data and other forms of 
communication. 

The Bell System has previously an- 
nounced plans to develop a satellite 
communication system and put it into 
operation as soon as possible. A 
ground station in Maine will be ready 
early next year. The first test satel- 
lite will be ready at the same time if 
the government designates a launching 
vehicle soon. The Bell System has 
offered to pay costs of the rocket 
and launching. 

Dr. Fisk said that the first experi- 
mental phone satellite will be roughly 
spherical, but may have as many as 
60 flat surfaces or gem-like facets. 
The solar cells on each facet would 
receive an equal amount of light. 

The solar cells—invented at Bell 
Laboratories—are small silicon wafers 
mounted on the outside surface to 
convert sunlight to electricity. The 
sapphire sheathing will protect the 
solar cells from the deteriorating ef- 
fects of electron bombardment. in 
space. It will also reduce the effects 
of proton bombardment. 

Two different sizes of test satellite 
—with diameters of 27 in. and 4 ft.— 
are being considered for experiments. 
More than half the surface of the 
smaller one would be covered by over 
4,000 Bell solar cells; the larger satel- 
lite, by 12,000 cells. 
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Covering the cells will be an equal 
number of slices of man-made, crystal- 
clear sapphire. Other protective ma- 
terials are not being ignored, although 
sapphire appears most promising. 

The sapphire does not cost signifi- 
cantly more than other high-quality 
materials that might be used, and 
only about one-third of the cost of 
the solar cells themselves. Sapphire 
has other advantages. Fer one, it 
will convey heat away, preventing the 
solar cells from being overheated in 
long periods of continuous sunlight. 

Another advantage is that the use 
of sapphire overcomes the effects of 
very rapid and extreme changes of 
temperature in space. Most materials 
expand when heated, but by different 
amounts. If the rate of expansion 
differs too much, the materials will 
not stay fastened together. 

The sapphire slices will be brazed 
to platinum sidebars which will be 
soldered to a ceramic base. These 
materials were selected because they 
all expand at very nearly the same 
rate. 

In addition, sapphire, being hard, 
will resist erosion by micrometeorites. 

Under the reliability approach 
evolved in the development of subma- 
rine cable and other systems, the goal is 
to assure a predictable long life for the 
entire system. To do this, engineers 
must develop each individual com- 
ponent to such a reliability that its 
individual predicted life-length is far 
greater than the goal for the over-all 
system. 

“First they must study the en- 
vironument—whether it is two miles 
beneath the sea or 4,000 miles above 
it—and define all the conditions and 
hazards,” Dr. Fisk said. 

Then, in an artificially simulated 
environment, they intensify the effects 
of the hazards to make the components 
fail. The engineers then redesign 
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until they can predict accurately the 
reliability of a required component. 
Finally, they conduct a real-life test 
program to confirm that the system 
performs as predicted. 

Unfortunately, the space environ- 
ment is not yet known precisely. The 
exact amounts of radiation in space 
at various altitudes and under varying 
circumstances have not been meas- 
ured. Therefore Bell Laboratories 
engineers plan to use their first ex- 
perimental satellite both to take ad- 
ditional measurements of space radia- 
tion and to serve as a first ‘‘real-life 
test” of components planned for a 
working satellite system. 


Heat-Shrinkable Tubing. The de- 
velopment of a radically new type of 
tubing that shrinks to form a tight 
mechanical bond within 7 seconds fol- 
lowing application of heat (135° C) 
has been announced by Alpha Wire 
Corp., New York, N. Y. 

Marketed under the name of Alph- 
lex Shrinkable Tubing, the new tubing 
is an irradiated, flame-retardant, 
thermally stable, modified polyolefin 
base product, with applications limited 
only to the imagination of the design 
engineer. It is manufactured in an 
expanded form so that it may be easily 
slipped over wire and cables, hosing, 
terminals or conduit. When _ heat 
of approximately 135 C. is applied, 
the material shrinks to a reduced, 
predetermined diameter forming a 
secure, tight-fitting bond.  Longi- 
tudinal shrinkage is less than 10 per 
cent. To supply the heat the use of a 
heat gun is recommended. However, 
satisfactory results may be obtained 
by oven heating, radiant heat, dipping 
in hot liquid, ovens, soldering iron or 
open flame. 

Alphlex Shrinkable Tubing is mar- 
keted in standard packages of 4-ft. 
lengths. Sizes 24 through 4 are avail- 
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able in black, white, red or yellow; 
sizes 2 through 1 (in.) available in 
black; other colors and lengths may 
be obtained on special order. 

When used as an insulating material, 
Alphlex Shrinkable Tubing has excel- 
lent temperature characteristics and 
remains flexible and strong from —55 
C. to 135 C. continuous. Within this 
temperature range it will not melt, 
harden, run, crack, or blister. The 
high flexibility of Alphlex Shrinkable 
Tubing remains constant throughout 
the expansion and reduction process. 
Adding to its application as an insu- 
lating material is the high dielectric 
and mechanical strength of Alphlex 
Shrinkable Tubing and its resistance 
to fungus growth and moisture absorp- 
When used in hydraulic pro- 
cesses, Alphlex Shrinkable Tubing 
provides a flexible, chemical-resistant 
jacket or coupling. 

Other applications for Alphlex 
Shrinkable Tubing include use as in- 
sulation for cables, harnesses, pigtails, 
bus bars and electronic components. 
Alphlex Shrinkable Tubing also may 
be used to protect terminals, to mend 
or join shafts, as a convenient insula- 
tion for tool handles, and for perma- 
nent, easy to use, cable markers and 
color coding.  Alphlex Shrinkable 
Tubing is easily marked or stamped 
with hot stamping equipment. 

Alphlex Shrinkable Tubing has an 
exceptionally long shelf-life and will 
retain its expanded diameter until 
heated. No radioactivity remains in 
or on the material after irradiation. 
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Private Satellite Telephone.— You 
could have a private telephone ‘‘line”’ 
between your own location and your 
business headquarters or friends’ 
homes many thousands of miles away 
with the space communication satel- 
lite system proposed by the General 
Electric Company. And the contact 
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could be made by a simple dialing of 
a telephone as you do now in the home 
or office. 

This feature of the proposed system 
was explained by A. G. Steinmayer, 
Jr., systems project engineer at the 
Company's Missile and Space Vehicle 
Department in Phildelphia, in a paper 
presented at: the National Symposium 
on Global Communication there. 

Mr. Steinmayer said an_ isolated 
business office in Australia could lease 
the service on a monthly basis and 
have a constant contact with London 
or with another office 75 miles away if 
desired. The satellite system would 
provide high quality reception at both 
terminals because the signals would 
not be affected by cloud cover or ad- 
verse weather conditions. 

A person living in an isolated area 
wishing to use the system would be 
equipped with broadcasting and _ re- 
ceiving equipment that would “tune 
in” on the radio frequencies assigned 
to him. By picking up his telephone, 
and dialing, he would be in contact 
with the satellite in his area. His 
signals would then be relayed by the 
satellite to a neighbor's ground station 
or to a commercial ground station 
which could handle his call through an 
adjoining ground communication 
system. 

The proposed General Electric sys- 
tem would employ ten or more satel- 
lites to cover the earth’s populated 
areas. Each satellite would have 
600 two-way voice channels. As one 
satellite passed out of range, another 
satellite would take over. Signals 
would go from a ground station to the 
satellite and down to a ground station. 
If the call were placed to the other 
side of the earth, the signals would go 
up to a satellite and back to a ground 
station several times. This enables 
the same channels to be used in at 
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least nine different areas of the earth 
and makes it possible for thousands of 
users to enjoy the benefits of satellite 
communication service. 

A major advantage of the General 
Electric system is that is could be 
used by all—small as well as large 


countries and private individuals. 
The cost for an isolated private indi- 
vidual would be high compared to an 
average urban phone bill but quite 
moderate, in terms of the service pro- 
vided, for an isolated user. 

As proposed by General Electric, 
its system is designed to benefit the 
entire populated world, rather than 
just the heavily populated and de- 
veloped areas. 

General Electric recently organized 
company —Communication 
to operate a space sat- 


a new 
Satellites, Inc. 
ellite communication system and has 
applied to the Federal Communica- 
tions Commission for an operating 
license. This organization would lease 
either individual channels or blocks 
of channels to users. 

The lease arrangement vastly simp- 
lifies the monitoring and billing prob- 
lem and completely obviates the need 
for the normal type of switchboard 
operation. 

In operation, each user would have 
a range of frequencies assigned to him 
that would carry the radio signals. 
Messages would be carried simultane- 
ously as they are on telephone wires 
today. This is technically known as 
frequency stacking messages. 

The system that is currently con- 
sidered to provide optimum benefits 
employs pulse code modulation tech- 
niques. This is transmission of the 
voice signals in pulses that are decoded 
at the receiving end. Full tonal 
fidelity is preserved so that the user’s 
voice loses no personality or inflection 


characteristics. 
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Machine Detects One Odd Atom in 
Crowd of a Billion.—A 4-ton research 
machine that detects one foreign atom 
lurking among more than a _ billion 
atoms—a feat comparable to picking 
out one wanted man hiding in a crowd 
of all the people in North and South 
America, Europe, and Africa—has 
been put to work at RCA’s David 
Sarnoff Research Center in Prince- 
ton, N. J. 

The new instrument, known as a 
spark-source mass spectrograph, is one 
of the world’s most sensitive research 
tools. It is being used by RCA scien- 
tists to probe new electronic materials 
that will be used in tomorrow’s com- 
puters, communications, and control 
systems. Built for RCA Laboratories 
by Associated Electrical Industries 
(Manchester), Ltd., in Manchester, 
England, the machine was delivered 
by transatlantic air freight for instal- 
lation at the Princeton center. 

Announcing the start of research 
operations with the new instrument, 
H. W. Leverenz, Director of Research, 
RCA Laboratories, said: 

“Basic research into the nature and 
properties of materials that have po- 
tentially useful electronic character- 
istics has become a principal founda- 
tion for the advances that we antici- 
pate for the future in electronics. 
RCA’s extensive and expanding pro- 
gram in materials research already 
has produced many results of both 
fundamental practical —im- 
portance. These materials make pos- 
sible new or improved devices now 
entering use in computers, television 
systems, microwave communications, 
satellites, and a multitude of defense 
electronic systems. This new pre- 
cision instrument, added to the ad- 
vanced facilities now in use at RCA 
Laboratories, will enable our scien- 
tists to probe even more deeply into 
the composition of promising new 
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materials, and thus to accelerate the 
development of new building blocks 
for tomorrow's electronic devices and 
equipment.” 

Mr. Leverenz described the 
mass spectrograph as an instrument 
whose sensitivity is ‘unmatched by 
that of any such equipment known to 
be in use in this country.”’ Describ- 
ing its ability to pick out and identify 
in only an hour, one impurity atom 
from among more than a billion atoms, 
he pointed out that even this minute 
degree of variation from a desired 
composition can affect the electronic 
properties of a material. 

The operating principle of the mass 
spectrograph is based on the fact that 
each chemical element has a distinct- 
ive mass and ionization characteristic 
distinguishing it from all other ele- 
ments. In the spark-source mass 
spectrograph, a small amount of the 
material to be analyzed is converted 
by an electric spark discharge into a 
cloud of positively-charged atoms, or 
positive ions. The charged atoms are 
accelerated, then sent in a narrow 
beam through electrical and magnetic 
fields which deflect them in varying 
degrees, depending upon their velocity 
and mass. The atoms of different 
elements are thus separated, to make 
impact at different points along a pre- 
cisely calibrated scale on a_photo- 
graphic plate. The result is a detailed 
record of all the elements present even 
in extremely minute quantities in the 
material that is being analyzed. 


new 


Publications Available from Library 
of Congress.—-The Library of Con- 
gress has a duplicate stock of back 
numbers of its publications, the East 
European Accessions Index and the 
Monthly Index of Russian Accessions, 
available for donation to libraries of 
educational institutions and public 
bodies in the United States. Serving 
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as keys to new literature coming into 
this country from the USSR and East 
Europe, these publications are of 
value to scientists and other scholars. 
They are more than mere ‘‘booklets,” 
for they contain information about 
the content of new Russian and East 
European books and articles. Each 
index gives in English a monthly ac- 
count of new material in a variety of 
subject fields as received from the 
USSR and East Europe by the Library 
of Congress and by other American 
research libraries. The translation of 
all books and articles into English 
and the elaborate subject guides in 
English permit a researcher who has 
no command of Russian or East 
European languages to identify easily 
the material important to him. 

The number of copies available of 
individual issues varies for both of 
these monthly publications, but the 
Library will undertake to fill requests 
for nearly complete back files insofar 
as the supply permits. All requests, 
whether for single issues or for all 
numbers available, will be filled in 
order of receipt. The requests should 
be addressed to the Chief, Exchange 
and Gift Division, Library of Cong- 
gress, Washington 25, D. C. 

Subscriptions for future issues of 
these publications may be placed with 
the Superintendent of Documents, 
United States Government Printing 
Office, Washington 25, D. C. 


“Plug Nozzle’? Rocket Engine.-A 
new type rocket engine whose design 
promises substantial cost and time 
savings in narrowing the gap between 
American and Soviet rocket propul- 
sion systems has been test fired by 
General Electric. 

Called the ‘‘plug nozzle” engine and 
termed the “first radical change in 
rocket engine design in more than a 
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decade,” the powerplant developed 
more than 50,000 Ib. of thrust in a 
recent firing at General Electric’s 
Rocket Test Station in Malta, N. Y. 

The company’s propulsion engineers 
claim the engine can be scaled up to 
multi-million-pound thrust levels in 
less time and at less cost than today’s 
conventional rocket propulsion 
systems. 

The feasibility of the G-E pioneered 
plug nozzle rocket engine concept, 
under study for several years at the 
firm’s Flight Propulsion Laboratory 
Department, is being demonstrated 
under contract with the National 
Aeronautics and Space Administra- 
tion. Recent firings of the engine 
are part of an investigation being 
conducted by the space agency. 

Dr. Kurt Berman, head of G.E.’s 
Rocket Test Station, said that the 
plug nozzle firing represents a major 
advance in the development of liquid 
propellant rocket engines of uncon- 
ventional design. plug nozzle,” 
he said, ‘‘is the first radical change in 
rocket engine design in more than a 
decade. We believe this engine can 
materially increase reliability and re- 
duce the time and cost of developing 
high-thrust rocket engines.” 

The plug nozzle engine differs from 
the conventional rocket engine both 
in appearance and development ap- 
proach. While existing rocket engines 
have an exhaust chamber which is 
bell-shaped, the plug nozzle is only 
half as long and is shaped like a funnel. 

According to G. E. engineers, in- 
creasing the thrust of conventional- 
type engines requires that the propel- 
lant injector be enlarged. However, 
this creates critical combustion prob- 
lems which can only be overcome with 
a major engineering and development 
effort. The results are costly, both in 
terms of long development cycles, and 
expenditures for constantly changing 
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hardware and _ high-thrust _ test 
facilities. 

The plug nozzle approach to higher 
thrust levels eliminates the need for 
“cut and try” methods of develop- 
ment. Combustion in the new engine 
occurs in small segmented combustor 
cells located around the periphery of 
the conical plug. These small cells 
effectively reduce the 
problem, and increase engine relia- 
bility. By increasing or decreasing 
the number of these cells, engines of 
various thrust levels can be built 
without developing a new injector for 
each engine size. 

“We can test segments of these 
cells,’ said Berman, “‘and then group 
the segments together to get as much 
thrust as we want. We don’t have 
to take the intermediate steps.” 

The plug nozzle is being studied for 
use as a multi-million-pound booster 
for large payloads, as an upper-stage 
engine and for space station propul- 
sion. It has a decided advantage over 
conventional engines in that varying 
the pressures in the combustor seg- 
ments affords a simple means to 
“steer” the rocket. 

Conventional engines, on the other 
hand, require a bulky gimballing 
structure to do the same thing. The 
gimballing structure mechanically 
turns the conventional engine to 
change the direction of thrust in 
respect to the direction of the rocket. 


combustion 


United Kingdom’s VERA Reactor. 
—The major purpose of the new 
VERA reactor at the Atomic Weapons 
Research Establishment, Alder- 
maston, England, which went critical 
for the first time late in February, is 
to gather information in an area where 
it has up to now been substantially 
lacking, that is, that of medium-en- 
richment, of diluted and of slightly 
moderated fast reactor systems. 
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6 ft. 


has walls 
thick and a ceiling 3 ft. thick. Around 
the cell are located the control room, 
the fuel element assembly room, labo- 
ratories, plant room for the ventila- 
tion system and filters and a labora- 
tory to house a 600-kV accelerator, 


The reactor cell 


still under construction. This ac- 
celerator will be used to produce deu- 
teron beams aimed at targets in sub- 
critical VERA cores which will yield 
bursts of neutrons in a series of kinetic 
experiments to measure reactivity in- 
crements and neutron lifetimes and to 
study neutron spectra with varying 
core compositions. 

Versatility is contributed to the 
reactor, and its name, (VERA— 
Versatile Experimental Reactor As- 
sembly) by the fuel element, which 
could hardly be more simple. Each 
element consists of two stainless steel 
U-channels into which }-in. thick, 
roughly square ‘‘biscuits’’ of enriched 
uranium, natural uranium and graph- 
ite are stacked in numbers and 
patterns which can be selected to cor- 
respond to practically any fuel dilu- 
tion or degree of moderation. 

Only the central portion of each 
element is occupied by this mixture, 
the remainder of the stack consisting 
of thicker blocks of natural uranium 
or other materials serving as the 
reflector. 

A small gap is left between the two 
channels when they are clamped to- 
gether so that the fuel element make- 
up can be checked at a glance—the 
enriched (93°) biscuits are coated 
with orange lacquer and the natural 
ones with green. This lacquer also 
helps to reduce the spread of radio- 
active dust. 

The reactor is brought into opera- 
tion by raising eight fuel elements, 
so that their central enriched por- 
tion enters the core. They are sus- 
pended by steel cables from winding 
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drums equipped with electromagnetic 
clutches. On quick release, all control 
and safety elements fall, out of the 
core, under gravity, closing the de- 
vice down immediately. 

Closed circuit TV is used to watch 
operations from the control room. 


Computer to Serve Advanced Tech- 
nologies.—A new scientific computer, 
designed to meet the advanced data 
processing needs of complex tech- 


nologies, has been announced by 
International Business Machines 
Corporation. 


The IBM 7072 will be able to proc- 
ess magnetic tapes recorded during 
engineering tests—jet aircraft wind 
tunnel trials, for example—despite 
a difference in the “‘languages’’ of the 
test cell and the computer. Binary 
data are entered directly into the com- 
puter, where they are converted into 
decimal form automatically and proc- 
essed for immediate use by engineers. 

The new system is suited especially 
for problems with an emphasis on 
powerful computing ability, rather 
than high-volume input and output 
of data. 

The latest member of a fully-com- 
patible family of computers, which 
includes the IBM 7070 and 7074 data- 
processing systems, the 7072 is ex- 
pected to serve a wide range of users, 
including aircraft manufacturers, en- 
gineering organizations and heavy 
industry. 

“At a time when science and indus- 
try are placing increased reliance on 
computing as they push into new 
fields, the 7072 offers great potential 
for intermediate-cost, high-powered 
technical data processing,” Y. P. 
Dawkins, director of marketing for 
IBM’s Data Processing Division, said. 
“The new system also provides broad 
business data processing capabilities.” 
Among the scientific and engineering 
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problems that will be handled by the 
7072 are: the creation of mathematical 
models to test product designs, sys- 
tems and operational theories ; simula- 
tion of space flight vehicle trajec- 
tories; design of nuclear reactors and 
supersonic aircraft; analysis of prod- 
ucts and structures to determine 
their ability to withstand stress; re- 
duction of large volumes of data to 
a form in which they can readily be 
analyzed; and solution of equations 
using such general mathematical tech- 
niques as matrix algebra and vector 
analysis. 

The 7072, with more than two times 
the internal speed of the IBM 704, 
has access to data stored in core 
storage memory in six-millionths of 
a second. 

The availability of 7070/7074 auto- 
matic programming systems will per- 
mit users to test their programs at 
any IBM Datacenter prior to the in- 
stallation of the 7072 and to put the 
system into immediate operation after 
it has been installed. These programs 
include IBM’s FORTRAN, AUTO- 
CODER, SORT and MERGE, RE- 
PORT GENERATOR and IOCS (In- 
put Output Control System). 


Device Analyzes Binary Gas Mix- 
tures.—The Gasalyzer, a new battery- 
operated portable instrument which 
accurately measures flue and auto 
exhaust gases, is now being marketed 
by the Industrial Equipment Division 
of Engelhard Industries, Inc., New- 
ark, N. J. 

The device, which weighs only 63 
lb., is capable of analyzing binary 
gas mixtures to within + 5 percent 
accuracy. Higher degrees of accuracy 
can obtained’ with special 
calibration. 

Developed originally for use by the 
U. S. Army Quartermaster Corps, 
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the Gasalyzer readily checks com- 
bustion efficiency and can be adapted 
to measure hydrogen, sulfur dioxide 
and carbon dioxide. It also can be 
used to measure fuel-air ratios in 
internal combustion engines, turbines 
and furnaces. 

Powered by a standard 6-volt, dry 
cell battery, the Gasalyzer will oper- 
ate as long as its battery has life. A 
dying battery will not affect perform- 
ance, and only when the electrical 
source is dead will the instrument stop 
functioning. During testing, bat- 
teries lasted for 72 continuous hours. 

The Gasalyzer can be adapted to 
110 volts a.c. by using a power pack 
available at additional cost. 

The operation of the Gasalyzer is 
based on the thermal conductivity 
principle using two sensory elements 
in a wheatstone bridge. The gas 
passes through an analysis cell, caus- 
ing the filament of the cell to heat or 
cool with respect to the filament of a 
reference cell. The temperature 
change of the analysis cell results in 
a resistance change which upsets the 
bridge balance. This change then is 
reflected on a highly sensitive meter. 


Photographic Exposure Repeater 
for Transistor Industry.—The first 
commercially available photographic 
exposure repeater able to produce ex- 
tremely precise photo-masks required 
by the transistor industry has been 
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developed by the David W. Mann 
Company, Lincoln, Mass., a division 
of the Geophysics Corporation of 
America. 

Built around a two-coordinate com- 
parator of the type which has identi- 
fied the company as a world leader 
in the very exact science of metrology, 
the new equipment permits position- 
ing of patterns to an accuracy of one 
micron (.0004 in.), exact focusing, 
and a precise ratio of reductions. The 
Repeater has a programming mecha- 
nism capable of performing various 
combinations of step-and-repeat op- 
erations for automatic exposures, al- 
though the instrument may also be 
operated manually for specially spaced 
exposures. 

A major feature of the new Re- 
peater is a three-pin plate positioning 
mount on the plateholder which per- 
mits accurate positioning ex- 
posure of an aligned series of plates. 
A mask use fixture of identical con- 
figuration makes it possible to insert 
a series of photo-masks and produce 
precisely aligned successive exposures 
of semiconductor material. 

The repeater is available in three 
basic models, an Industrial, Research, 
or Manual model. The instrument’s 
major advantages lie in cost, time 
savings, and proprietary security in 
the production of masks, and con- 
venience in experimentation with new 
types of masks. 
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This engineer is measuring the earth’s veloc- 
ity with respect to the sun. He is using the 
method of template spectroscopy to measure 
Doppler spectral shift in the optical range. 


This shift, which is due to the familiar Doppler effect that relates fre- 
quency and velocity, indicates a change in the relative velocity between the 
observing point on earth and the observed point on the sun. Studies have 
thus far demonstrated that with such an optical Doppler measuring system, 
we are able to measure velocity changes with a sensitivity of the order of 
meters per second. 


This is a portion of a study being made at The Franklin Institute 
Laboratories for the Air Force, whose purpose it is to determine whether 
natural radiation from the sun, stars and interstellar space might be used in 
space navigation. The work mentioned above is now being extended to 
consider stellar sources whose light energy is considerably less than that 
of the sun, and for which special detection techniques are required to 
derive useful signals. 


Velocity measurement by way of optical Doppler techniques is one of a 
variety of interesting research programs now being carried out for industry 
and government agencies in our Laboratories. To obtain additional in- 
formation concerning activities or staff positions in the Laboratories, address 
Director of Personnel, The Franklin Institute. 


THE FRANKLIN INSTITUTE LABORATORIES 
20th and Parkway, Philadelphia 3, Pennsylvania 


Chemistry—Physics—Solid State Sciences—Electrical Engineering—Mechanical and Nuclear 
Engineering—Engineering Psychology—Operations Analysis 
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Now in the fourth year of operation, The Franklin Institute Com- 
puting Center has shared its technical know-how, modern electronic 
equipment and experienced personnel with hundreds of progressive 
industries and government agencies across the nation. 


We offer the services of creative people, highly skilled in their respec- 
tive fields and ably trained in the application of these skills to the ever 
expanding area of electronic computers and data processing systems. 


This staff is now available for analysis, system design, programming 
or coding of projects of unlimited scope or context. Input to our large 
scale computer and peripheral equipment is acceptable in any form. 
Results are provided on cards, plastic or metallic tape, and in com- 
pletely edited printed copy. Our extensive library of automatic coding, 
engineering, data processing and mathematical routines is available to 
all users, and machine time is provided with or without the services 
of programming personnel. 


THE FRANKLIN INSTITUTE 
Computing Center 


20th Street and Benjamin Franklin Parkway 
Philadelphia 3, Pa. LOcust 4-3600 
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